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Time  to  go  before  release  command 
INS  gimbal  angle, 

State  transition  propogatlon  matrix,  and 
Roll  attitude  angle 
State  transition  matrix 
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LIST  OF  SYMBOLS  —  (continued) 


SYMBOL 


0 

0 


12 

13 


0) 


CG 


0) 


(1) 


tm 

sim 


DEFINITION 

State  transition  matrix 
State  transition  matrix 
Angular  rate 
Attacker  angular  rate 
Predicted  angular  rate,  and 
LOS  rate  component 

Commanded  angular  rate  for  gun 
LOS  angular  rate  component 
Angular  rate  of  glmbal  axis 
Angular  rate  of  target  LOS ,  and 
Component  of  LOS  rate 

Angular  rate  of  LOS  about  £  axis 
Angular  rate  of  LOS  about  m  axis 
Angular  rate  of  tracker  coordinates 
Component  of  angular  rate  vector,  w 
Component  of  angular  rate  vector,  w 
Component  of  angular  rate  vector,  u) 
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LIST  OF  SYMBOLS  -  /continued) 


SYMBOL  DEFINITION 


The  coordlnat'a  systems  used  are  right-handed  and  consist  of  the  following: 


N,  E,  D 


u,  V,  w 


u  t 
w 


w 

w 


B,  Z,  m 


t,  c,  d 
®h’  “h 


Inertial  or  earth-fixed  coordinates; 

N  -  north,  E  -  east,  and  D  -  down. 

Attacker  body  coordlntes; 

u  -  forward,  v  -  out  to  right  wing,  and  w  -  down. 
Attacker  wind  coordinates; 

u  -  along  attacker  velocity  relative  to  airmass, 

V  -  to  the  right,  and  w  -  down. 

Sight  coordinates; 

8  -  along  the  target  line-of-slght, 

Z  -  out  to  the  right,  and  m  -  down. 

Tracker  coordinates,  t  -  along  the  tracker  axis, 
c  -  to  the  right,  and  d  -  down. 

Helmet  sight  coordinates,  Sj^  -  along  the  target  sight, 
out  to  the  right,  and  m^  -  down. 


Superscripts 

Estimate 

.  Derivative 

_  Vector  quanity 

.p  Vector  transpose 


1.0  INTRODUCTION 


The  Advanced  FIREFLY  Assessment  (AFFA)  Program  extends  and  generalizes 
the  Integrated  Fire  and  Flight  Control  (IFFC)  studies  so  effectively  performed  under 
the  FIREFLY  II  Program.  These  studies  have  shown  the  potential  for  significant 
Increases  In  aircraft  survivability  and  effectiveness  offered  by  Integration  of  the  fire 
control  and  flight  control  functions  during  the  process  of  delivering  a  weapon  to  Its  tar¬ 
get.  To  date,  these  studies  have  provided  Initial  concepts  and  preliminary  systems 
definitions.  The  AFFA  Program  will  extend  this  woric  to  operational  applications,  with 
specific  avionics,  and  therefore  requires  a  more  versatile  analysis  tool  for  In-depth 
analysis. 

The  Advanced  FIREFLY  Assessment  Program  has  as  Its  primary  goal  the  design 
of  a  Generalized  Mechanization  of  FIREFLY  concepts,  the  development  of  a  digital  sim¬ 
ulation  of  the  application  of  this  Generalized  Mechanization  to  a  given  aircraft  (e.  g. , 
F-16)  operating  in  a  given  mode  (e.g. ,  air-to-air  gunnery),  and  the  evaluation  of  the 
resulting  designs.  Generalization  as  used  here  Implies  the  ability  of  each  function  In 
the  mechanization  to  encompass  different  avionic  subsystems  operating  In  different 
aircraft  In  different  modes.  This  Generalized  Mechanization  will  also  allow  the  effec¬ 
tive  Investigation  of  Advanced  FIREFLY  applications  beyond  those  cases  considered  In 
the  FIREFLY  II  Investigations. 

Although  the  primary  thrust  of  the  AFFA  Program  involves  determining  a  pre¬ 
ferred  fire  control  technology,  associated  avionics  subsystem  requirements  and  con¬ 
trol  law  mechanizations,  it  Is  not  possible  to  achieve  the  most  effective  Integration  with 
the  flight  control  system  without  consideration  of  aircraft  and  flight  control  character¬ 
istics,  Accordingly,  the  AFFA  Program  shall  establish  the  extent  to  which  consider¬ 
ation  of  detailed  aircraft/flight  control  dynamics  must  be  implemented,  short  of  control 
system  redesign.  The  goal  here  Is  clearly  overall  system  performance  as  this  is 
determined  by  dynamic  Interaction  of  the  major  functions  and/or  subsystems. 

The  AFFA  Program  is  particularly  timely  because  continuing  r^ld  advances  in 
microelectronics  and  the  imminence  of  both  Very  Large  Scale  Integration  (VLSI)  and 
Very  High  ^eed  Integration  (VHSI)  make  very  large  Increases  In  computer  capability 
both  feasible  and  affordable.  Wlih  such  increased  processing  and  memory  capacity, 
much  more  sophisticated  and  complex  algorithms  can  be  mechanized  to  broaden  the 
vistas  of  operational  optimizations.  As  a  result,  the  FIREFLY  aircraft  will  be  capable 
of  more  effective  maneuvering  prior  to  as  well  as  during  weapon  release. 
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Currently  each  aircraft  with  its  avionic  sensors  and  associated  flight  control  and 
fire  control  logic  must  be  separately  analyzed,  simulated,  and  evaluated  for  its  per¬ 
formance  in  a  given  tactical  situation.  When  the  mission  scenario,  aircraft,  or  sensor 
combination  changes,  extensive  and  costly  effort  results  in  arriving  at  a  workable 
flight/fire  control  mechanization  that  is  properly  designed  and  evaluated.  This  diffi¬ 
culty  motivates  a  generalized  and  modularized  characterization  of  an  integrated  flight/ 
fire  control  system  mechanization  and  associated  algorithms.  The  idea  of  such  a 
Generalized  Mechanization  is  to  Incorporate  common  functions  (l.e.,  atmospheric  and 
wind  models,  target  models,  attacker  models,  estimation  algorithms,  control  logics, 
etc.)  so  as  to  not  have  to  start  from  scratch  each  time  the  aircraft,  sensor,  and/or 
mission  scenario  change. 

The  term  "Generalized  Mechanization"  will  mean  sets  of  equations,  algorithms, 
and  control  laws  capable  of  being  programmed  into  an  airborne  digital  computer  which 
can  be  specialized  to  perform  successfully  in  a  variety  of  tactical  aircraft  with  differ¬ 
ing  avionics  sensors,  fire  control  algorithms,  control  laws,  flight  control  systems, 
and  weapon  delivery  modes. 

The  design  of  the  Generalized  Fire  Control  Mechanization  is  of  critical  impor¬ 
tance  to  the  success  of  the  AFFA  Study.  If  the  Generalized  Mechanization  is  not  suf¬ 
ficiently  flexible,  considerable  duplication  of  effort  may  result  during  the  development 
of  the  FIREFLY  Simulation  and  Analysis  Program  and  also  during  the  design  effort 
called  for  in  Task  4.  The  Generalized  Mechanization  permits  a  unified  approach  to  be 
taken  in  the  overall  avionics  integration  without  getting  umiecessarily  bogged  down  in 
technical  detail.  Looking  beyond  the  immediate  alms  of  the  FIREFLY  Assessment 
Study,  the  generalized  nature  of  the  software  architecture  will  be  of  great  value  in 
implementing  tlie  stale  estimator  and  fire  control  algorithms  in  an  actual  operational 
flight  program.  Therefore,  a  great  deal  of  emphasis  has  been  placed  in  developing  a 
set  of  requirements  for  the  fire  control  mechanization  that  has  maximum  generality 
with  regard  to  sensor  inputs  and  weapon  delivery  modes. 

This  report  documents  the  requirements  to  which  the  AFFA  Generalized 
Mechanization  is  to  be  designed.  This  is  done  by  first  addressing  the  scope  of  the 
Requirements  Analysis  followed  by  a  description  of  the  analysis  done  in  deriving  the 
requirements  of  the  generalized  mechanization  in  Section  3.  Finally,  the  requirements 
are  summarized  in  Section  4  discussing  them  explicitly  in  relation  to  their  place  in  the 
overall  fire  control  system. 
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A  discussion  of  die  work  done  on  the  FIREFLY  II  familiarization  and  Initial 
analysis  of  two  of  the  potential  advanced  concepts  for  evaluation  is  contained  in  the 
appendices. 


2. 0  GENERALIZED  MECHANIZATION  REQUIREMENTS  ANALYSIS  SCOPE 


2.1  GOALS  AND  SCOPE  OF  THE  GMRR 

The  goals  and  scope  of  the  Generalized  Mechanization  Requirements  Report 
(GMRR)  are  to  specify  detailed  and  explicit  requirements  for  the  Generalized  Mechan¬ 
ization  portion  of  a  software  tool  called  the  FIREFLY  Simulation  and  Analysis  Program 
(FSAP).  In  this  regard  the  GMRR  is  restricting  the  statement  of  GM  requirements  to; 

1.  An  identification  of  a  set  of  functions  which  comprise  the  Generalized  Mech¬ 
anization.  This  essentially  characterizes  the  gross  structure  of  the  Gener¬ 
alized  Mechanization  and  identifies  specific  Interfaces. 

2.  An  identification  of  all  inputs  and  outputs  of  each  function  that  is  sufficient 
to  encompass  the  major  weaponry  modes  of  air-air  gunnery,  air-ground 
gunnery,  and  bombing. 

3.  A  specification  of  all  interface  requirements  between  functions  and  con¬ 
straints  on  ir^uts /outputs  of  each  function  of  the  Generalized  Mechanization. 

4.  An  identification  of  all  algorithms  and  their  associated  functional  require¬ 
ments  for  each  Generalized  Mechanization  function. 

5.  An  identification  of  coordinate  frame  options  and  state  variables  for  each  of 
the  three  weaponry  modes. 

6.  An  outline  or  preliminary  statement  of  the  mathematical  equations  which 
characterize  the  Generalized  Mechanization. 

A  detailed  specification  of  the  error  model  requirements  will  be  included  as  part 
of  the  Avionic  Sensor  functional  requirements.  A  key  issue  v/hich  limits  the  scope  of 
the  GMRR  is  the  boundary  between  the  GMRR  and  the  Generalized  Mechanization  Design 
Report  (GMDR).  In  order  to  clarify  this  boundary,  the  scope  of  the  content  of  the 
GMDR  is  presented  below  thereby  showing  how  the  combination  of  the  GMRR  and  the 
GMDR  feed  into  structuring  the  FIREFLY  Simulation  and  Analysis  Program. 

GMDR  SCOPE 

1.  Evaluation  and  choice  of  coordinate  frrmes  for  different  avionic  subsystem 
combinations  and  different  weaponry  /nodes. 

2.  Evaluation  and  choice  of  state  variables  for  different  weaponry  modes. 


4 


3.  Evaluation  and  choice  of  algorithms  within  each  function  of  the  GM.  This 
includes  detailed  statement  of  all  equations  to  be  implemented  in  the 
FIREFLY  Simulation  and  Analysis  program. 

4.  Identification  of  how  the  inputs  are  determined  from  the  inputs  for  each 
function. 

5.  Justification  of  each  assumption  and  identification  of  the  limits  on  various 
parameters  and  variables  to  yield  specific  requirement  within  each  GM 
function. 

6.  Identification  of  all  assumptions  made  for  each  algorithm  of  each  function. 

The  Generalized  Mechanization  (GM)  consists  of  sets  of  equations,  algorithms, 
and  control  laws  capable  of  being  programmed  into  an  airborne  digital  computer  which 
can  be  specialized  to  perform  successfully  in  a  variety  of  tactical  aircraft  with  differ¬ 
ing  avionics  sensors.  Figure  1  illustrates  how  the  GM  relates  to  the  aircraft  and  to 
the  other  elements  of  the  tactical  environment. 

Requirements  on  the  GM  will  be  Imposed  by  the  following  considerations: 

1.  the  mathematical  adequacy  of  the  system  dynamics  formulation, 

2.  software  requirements  and, 

3.  hardware  requirements. 

Each  of  these  will  result  in  constraints  and  limitations  on  the  GM.  The  current 
version  of  this  document  focuses  on  the  mathematical  adequacy  of  the  GM  and  on  the 
variety  of  sensors  addressed  by  the  GM.  Requirements  imposed  upon  the  GM  due  to 
the  other  considerations  will  be  determined. 

Figure  2  illustrates  the  development  flow  for  the  Generalized  Mechanization 
Requirements.  Section  3. 0  of  this  report  wlU  discuss  the  Generalized  Mechanization 
Requirements  Analysis  and  Section  4. 0  will  contain  a  complete  presentation  of  the 
resulting  requirements. 


TARGET  pilot 


FIGURE  1.  FUNCTIONAL  FLOW 
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FIGURE  2.  REQUIREMENTS  DEVELOPMENT  FLOW 


3.  0  GENERALIZED  MECHANIZATION  REQUIREMENTS  ANALYSIS 


The  Generalized  FIREFLY  Mechanization  is  to  be  directly  applicable  to  three 
basic  weapon  delivery  modes;  air-to-air  gunnery  (AAG),  air-to-ground  giuuiery  (AGG), 
aiid  bombing.  In  addition,  it  is  to  have  inherent  compatibility  with  more  conventional 
weapon  delivery  modes  such  as  blind  coordinate  bombing,  offset  bombing,  angle  rate 
bombing,  and  all  other  types  of  Continuously  Computed  Release  Point  (CCRP)  weapon 
delivery  modes. 

The  AAG  and  AGG  inodes  are  director  gunfire  control  modes  in  which  the  predic¬ 
tion  of  future  target  motion  is  made  on  the  basis  of  data  from  a  tracking  sensor,  such 
as  a  radar  or  electro-optical  tracker.  Based  on  range,  range  rate,  angle,  and  angle 
rate  measurements  from  the  sensor  or  sensors  combined  with  ownship  data,  estimates 
of  target  position,  velocity  and  acceleration  are  formed  in  the  estimator  fimction  of 
the  ^’^neralized  Mechanization  (GM).  These  target  state  estimates  are  used  in  the  fire 
control  function  of  the  GM  to  form  an  extrapolated  estimate  of  the  target's  position  at 
some  future  time.  Simviltaneously,  the  fire  control  function  computes  the  ballistic 
trajectory  of  the  projectile  out  to  one  time-of-flight  based  on  projectile  aerodynamic 
forces  and  gravity. 

The  bombing  mode  is  based  on  CCRP  after  target  designation.  Target  designation 
can  be  accomplished  either  by  lock-on  with  a  tracking  sensor  or  by  overlaying  a  dis¬ 
played  target  with  a  designation  symbol  and  "marking"  its  position  at  tuo  instant  of 
designation.  The  target  position  vector  after  designation  is  computed  either  from 
tracker  inputs  or  from  aircraft  velocity  data,  hi  general,  this  velocity  data  may  be 
either  air  data  or  inertial  velocities.  The  impact  point  vector  and  the  predicted  future 
target  position  vector  are  computed  continuously  in  the  CCRP  modes.  When  the  differ¬ 
ence  between  these  two  vectors  is  reduced  to  an,  acceptable  level,  the  selected  weapon 
is  released.  The  fire  control  system  generates  the  control  laws  for  fuselage  aiming 
in  the  AAG  and  AGG  modes  and  the  steering  command  inputs  for  guiding  the  aircraft 
automatically  to  the  weapon  release  point  in  a  bombing  mission.  These  control  laws 
are  applied  to  the  flight  control  system  through  an  interface  (coupler)  to  achieve  the 
desired  aircraft  response  characteristics  to  command  angular  rates. 

3.1  FUNCTIONAL  RELATIONSHIPS  AND  SIGNAI;  FLOW  DIAGRAMS 

The  overall  signal  flow  diagram  of  the  geneialized  FIREFLY  control  system 
mechanization  is  shown  in  Figure  3  The  GM,  shown  in  dotted  lines,  receives  input 
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FIGURE  3,  SIGNAL  FLOW  DIAGRAM  OF  FIREFLY  IMPLEMENTATION 


from  the  avionics  subsystems  (sensors)  through  an  interface  and  generates  the  display 
information,  aircraft  steering  control  laws,  and  weapon  release  commands  for  various 
weapon  delivery  modes.  The  control  laws  in  turn  are  applied  through  an  interface  (or 
coupler)  to  the  aircraft  flight  control  system,  the  display  information  is  applied  through 
an  interface  to  the  A/C  pilot  interface,  and  the  weapon  release  commands  are  applied 
to  the  weapons  through  a  weapons  interface.  These  interfeces  perform  the  operations 
or  conversions  required  to  ensure  compatibility  between  the  FIREFLY  GM,  the  sen¬ 
sors,  and  the  flight  control  system.  These  interfaces  are  contained  within  an  Opera¬ 
tional  Flight  Program  (OFP)  and  are  external  to  the  GM. 

Figures  4  through  9  show  the  input /output  relationship  for  various  subfunc¬ 
tions  in  the  GM  and  the  interdependence  between  the  various  computational  units.  The 
avionic  subsystem  outputs  are  processed  through  state  estimators  which  compute  the 
estimated  quantities  needed  for  the  fire  control  solution.  The  ownship  estimator 
shown  in  Figure  3-2  combines  the  body  sensor  outputs  to  generate  smoothed  estimates 
of  aircraft  position,  inertial  velocity  and  acceleration,  angle  of  attack,  and  sideslip. 

The  atmospheric  estimator  (Figure  5  combines  the  air  data  measurements  with  the 
ownship  estimator  outputs  to  obtain  wind,  aircraft  velocity  relative  to  air  mass,  air 
density,  and  altitude  information.  These  estimators  may  range  in  complexity  from 
simple  mixing  devices  to  more  sophisticated  Kalman  filters  depending  upon  the  accur¬ 
acy  requirements  in  a  specific  application.  Ihe  target  state  estimator  shown  in  Fig¬ 
ure  6  uses  Kalman  filtering  algorithms  to  obtain  accurate  estimates  of  target  range, 
relative  velocity,  and  absolute  target  acceleration  from  ownship  state  estimates  and 
line-of-sight  (LOS)  sensors. 

Figures  7  and  8  show  the  functional  relationships  for  the  gxmnery  and  bomb¬ 
ing  modes.  Tlie  fire  control  solution  generates  the  gim  pointing  commands,  predicted 
future  target  position,  miss  distance,  and  other  variables  needed  for  the  gunnery 
modes.  The  time-of-flight  and  the  vectors  Rg  and  G  of  the  projectile  ballistic  which 
are  also  needed  for  the  fire  control  solution  are  obtained  from  the  ballistic  algorithms. 

The  functional  relationships  for  generating  the  bombing  solution  and  the  origin  of 
various  variables,  such  as  measured,  estimated,  or  known,  are  summarized  in 
Table  1.  The  variables  labelled  as  "known”  are  either  fixed  for  a  given  aircraft  or 
specified  by  the  pilot  to  achieve  the  desired  bombing  accuracy. 

The  control  laws  and  the  flight  control  system  interface  is  shown  in  Figure  9 
for  all  three  weapon  delivery  modes.  There  are  basically  two  control  laws  that  the 
GM  must  generate;  one  is  the  fuselage  aiming  commands  for  AAG  and  AGG  and  the 
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FIGURE  4.  INPUT/OUTPUT  REQUIREMENTS  FOR  THE  OWNSHIP  ESTIMATOR 
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FIGURE  6.  IKPUT/OUTPUT  REQUIREMSNTS  FOR  THE  TARGET  STATE  ESTIMATOR 
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FIGURE  8.  FUNCTIONAL  REQUIREMENTS  FOR  THE  BOMBING  FIRE  CONTROL  SOLUTION 
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FIGURE  9.  CONTROL  LAW  -  FLIGHT  CONTROL  SYSTEM  INTERFACE 


TABLE  1.  INPUTS  PJEQUIRED  FOR  A  BOMBING  SOLUTION 


SY^fflOL 

DESCRIPTION 

ORIGIN 

"c 

•k 

Ownshlp  inertial  velocity 

Ovmshlp  state  estimator 

"c 

Ownshlp  Inertial  acceleration 

Ownship  state  estimator 

"CA 

Ownshlp  velocity  relative  to  air 

Atmosphere  estimator 

w 

Wind  velocity  relative  to 
inertial  coordinates 

Atmosphere  estimator 

p 

Air  density 

Atmosphere  estimator 

Ve 

Bomb  ejection  velocity 

Known 

Sight  parallax  relative  to  body 

h 

Bomb  (stores)  parallax  relative 
to  body 

Known 

°of 

Target  position  offset  for 
stick  or  offset  bombing 

Known 

T 

Rack  delay 

Known 

\ 

Relative  target  velocity 

Target  state  estimator 

Inertial  target  acceleration 

II  II  It 

R 

Range  to  target 

II  It  It 

Azimuth  and  elevation  angles 
of  the  LOS  relative  to  body 
coordinates 

11  II  It 

u 

Ownship  angular  rate  vector 

Avionic  Subsystem  (sensors) 

[E] 

Transformation  matrix  from 
inertial  to  body  axes 

II  II  II 

h 

c 

Altitude  above  target 

"  "  (altimeter) 

Gravity  drop  along  the  local 
velocity  vector  relative  to  air. 

Ballistic  algorithm 

G' 

Ballistic  vector  along  the  local 
vertical 

II  It 

Time  of  flight 

II  It 

DGO 

Distance  to  go  before  turning 
(for  stick  bombing) 

Computed 

Time  to  go  before  bomb  release 

Computed 

(*)  Inertial  quantities  are  defined  in  Section  3.3.1 
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other  one  is  the  aircraft  steering  commands  for  bombing.  Additional  feedbacks  that 
are  required  to  essentially  decouple  the  lateral  and  longitudinal  responses  of  the  air¬ 
craft  for  the  AGG  mode  are  also  shown. 

The  GM  basically  performs  four  functions;  measurement  function,  state  estima¬ 
tion  function,  fire  control  solution,  and  the  generation  of  steering  command  signals 
(control  laws).  From  a  requirement  dependency  standpoint  the  flow  is  reversed  as 
shown  in  Figure  10.  Design  considerations  dictate  the  requirements  for  the  fire  con¬ 
trol  solution.  From  these  requirements  the  state  estimation  and  control  law  require¬ 
ments  are  determined.  Specific  sensors  on  board  an  aircraft  are  adapted  to  the 
generalized  state  estimator  and  to  the  fire  control  solution  thro’igh  the  measurement 
function.  Similarly,  the  control  law  requirements  are  met  on  a  specific  aircraft  (and 
a  given  flight  control  system)  by  proper  design  of  a  specialized  coupler. 

The  analysis  of  the  requirements  for  each  of  the  four  basic  functions  will  be  con¬ 
tained  in  the  following  sections  based  on  the  GM  requirements  for  AAG,  AGG,  and 
bombing  delivery  modes. 


FIGURE  10.  FUNCTIONAL  DEPENDENCY  FLOW  OF  THE  GENERALIZED  MECHANIZATION 


3.2  HRE  CONTROL  SOLUTIONS 

Several  refinements  and  improvements  for  a  generalized  fire  control  solution 
were  discovered  during  the  HREFLY  familiarization  phase.  There  refinements  are 
imposed  as  requirements  on  the  generalized  FIREFLY  mechanization  design,  sum¬ 
marized  in  Section  4. 

This  section  includes  a  discussion  on  ballistic  algorithms  (Section  3.2. 1),  time 
of  fli^t  computation  (Sections  3. 2. 3  and  3. 2. 5),  Fire  control  solutions  for  AA.G  and 
AGG  modes  (Section  3.2.2),  for  bombing  (Section  3.2. 1)  and  for  stick  bombing 
(Section  3.2. 6). 

3.2.1  Requirements  Analysis  For  Ballistic  Algorithms 

This  section  deals  with  the  definition  of  the  requirements  of  air-to-air  (A/A)  and 
air-to-ground  (A/G)  ballistic  algorithms  and  the  potential  utilization  of  some  of  these 
in  FIREFLY.  Section  3. 2. 1. 1  discvisses  the  general  requirements  for  A/G  weapon 
delivery,  followed  by  a  comparison  of  the  relative  merits  of  closed  form  ejq)ressions 
vs.  numerical  integration  ballistic  formulations.  The  requirements  for  A/A  fire  con¬ 
trol  are  discussed  in  Section  3.2. 1.2. 

3. 2. 1. 1  A/G  Ballistic  Algorithms 

For  the  purpose  of  bombing,  AGG  and  rocketry,  there  are  certain  require¬ 
ments  on  the  aircraft  in  order  to  accomplish  effective  weapon  delivery.  The  aircraft 
sensors  and/or  estimators  must  be  able  to  measure  or  estimate  the  relative  target 
position  and  velocitj ,  die  aircraft  velocity  in  the  air  mass,  and  the  gravity  vector. 

In  addition,  in  order  to  predict  where  the  weapon  will  fall,  there  exists  a  requirement 
to  have  available  in  the  airborne  fire  control  computer  a  ballistic  algorithm  for  the 
following  weapon  types:  low  drag  bombs,  retarded  bombs,  cluster  weapons,  gun  pro¬ 
jectiles,  and  unguided  rockets.  These  weapons  have  in  common  the  characteristics 
that  they  are  all  ballistic  projectiles,  that  is,  that  once  released  from  the  aircraft, 
the  only  forces  which  act  on  the  weapon  are  the  aerodynamic  forces  and  the  force  due 
to  gravity.  Aerodynamic  forces  usually  include  drag,  velocity  jump  and  windage  jump 
efforts.  (Note  that  for  the  case  of  unguided  rockets,  the  weapon  can  be  considered  a 
ballistic  projectile  only  after  burnout). 

In  the  past,  attempts  were  made  to  provide  a  ballistic  algorithm  in  closed- form  by 
formulating  a  simplified  drag  model  in  such  a  way  that  the  differential  equation  of 
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motion  is  solvable  by  analytical  methods.  In  addition,  these  approximate  solutions 
may  be  modified  by  empirically-derived  functions  that  are  Intended  to  compensate 
for  some  of  the  simplifications  in  the  drag  model,  order  to  provide  reasonable  ac¬ 
curacy  over  a  specified  release  envelope. 

The  other  method  of  providing  an  airborne  ballistic  algorithm  is  that  of  numerical 
integration.  While  this  method  is  advantageous  from  the  point  of  view  of  being  very 
accurate  and  flexible,  it  is  not  directly  amenable  to  the  type  of  fire  control  solution 
sou^t  in  the  ilREFLY  algorithm.  Adaptation  of  numerical  integration  methods  would 
require  an  iterative  solution  which  may  take  large  amounts  of  computer  time  (up  to  a 
few  seconds)  to  complete  the  calculations.  However,  it  is  possible  to  modify  tbe  num¬ 
erical  integration  technique  for  near-real-time  solution  in  an  airborne  computer.  The 
premimn  example  of  the  munerical  integration  technique  of  computing  the  weapon 
trajectory  in  an  airborne  ballistic  algorithm  is  the  China  Lake  Algorithm. 

The  following  is  a  brief  discussion  of  the  relative  advantages  of  the  closed-form 
and  numerical  integration  ballistic  algorithm  techniques. 

1.  Both  types  of  algorithms  are  nearly  equal  in  accuracy  when  a  single  weapon 
is  considered  in  the  "heart  of  the  delivery  envelope.  "  However,  imder  more 
severe  delivery  conditions  the  numerical  integration  technique  is  superior 

to  the  closed-form  method,  due  to  its  generality,  because  the  empirical 
modifying  functions  employed  in  the  closed-form  method  are  derived  only  for 
a  particular  region  of  the  delivery  envelope. 

2.  The  numerical  integration  method  is  more  flexible  vhea  a  new  weapon  is  to 
be  added  to  the  algorithm;  all  that  need  be  done  is  to  fit  the  aerodynamic 
functions  of  that  weapon  into  the  format  of  the  algorithm.  On  the  other  hand, 
for  the  closed-form  method  an  aitirely  new  function  must  be  derived  for 
eadi  new  weapon,  thus  leading  to  larger  memory  requirements. 

3.  The  one  primary  advantage  which  tlie  closed-form  method  has  over  the  num¬ 
erical  integration  method  is  that  of  computational  time.  While  the  munerical 
integration  method  is  more  accurate  over  a  larger  delivery  envelope,  and  is 
a  more  flexible  algorithm  than  the  closed-form  method,  the  larger  number 
of  calcula.tions  required  for  the  numerical  integration  makes  the  total  compu¬ 
tational  time  greater  than  that  for  the  closed-form  method. 

In  conclusion,  it  is  apparent  that  both  methods  have  their  own  merits,  and  that  when 
when  only  a  very  few  weapons  are  included  in  the  algorithm  the  two  methods  are  approx- 
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imately  equal  in  value.  However,  when  there  are  many  weapons  of  different  types  in 
the  aircraft  weapon  repertoire,  the  numerical  integration  method  is  most  likely  the 
more  effective  ol’  the  two  because  of  its  greater  flexibility,  accuracy,  and  the  less 
memory  space  ce-quired.  The  cost  of  memory  space,  however,  is  unlikely  to  be  a 
determining  factor  since  smaller  and  faster  computers  are  continuously  being 
developed. 

3. 2. 1. 2  A/A  Ballistic  Algoritiuns 

Just  as  in  the  case  of  A/G  weapon  delivery,  there  exists  a  requirement  in  AAG 
to  have  a  ballistic  algorithm  available  in  tiie  fire  control  computer.  Based  iq>oa  the 
knowledge  and  e:}q)erience  gained  with  A/G  ballistic  algorithms,  it  is  evident  that  an 
A/A  ninnerlcal  integration  ballistic  algorithm  would  probably  yield  greater  accuracy 
than  a  fitted  ballistic  algorithm  in  closed-form.  At  present  there  are  no  A/A  ballistic 
algorithms  which  make  use  of  the  numerical  integration  technique  similar  to  that  of 
the  China  Lake  Algorithm;  there  are  numerous  closed-form  ballistic  formulations  for 
AAG,  each  of  which  makes  different  simplifying  assumptions  according  to  which  gun- 
si^t  lead  ang^e  computation  the  ballistic  eqxiatims  require.  The  disadvantages  pre¬ 
sent  in  A/G  closed-form  ballistic  algorithms  —  namely,  the  large  memory  require¬ 
ments  brought  about  by  a  large  number  of  weapon  types,  and  the  inaccuracies  which 
occur  when  releases  are  made  outside  of  the  ’heart  of  the  delivery  envelcqpe"  —  az’e 
somewhat  alleviated  in  the  A/A  case.  Memory  requirements  would  be  smaller  be¬ 
cause  the  number  of  different  weapons  would  be  small  for  A/A  applications.  Also, 
in  A/A  encounters  the  ranges  involved  are  not  as  large,  and  the  analog  of  the  "release 
envelope"  is  not  as  severely  large  as  in  A/G  weapon  delivery.  All  of  this  puts  the 
closed-form  ballistic  algorithms  on  a  par  with  the  niunerical  integration  technique. 

In  summary,  even  though  the  numerical  integration  and  closed-form  ballistic 
algorithms  are  of  aH>roximately  the  same  value,  from  the  point  of  view  of  c'  oomon- 
ality  with  the  A/G  ballistic  algorithm,  the  numerical  integration  method  is  most 
likely  the  better  of  the  two  methods.  However,  caution  must  be  exercised  so  that 
t:q)date  rate  requirements  imposed  on  the  baUistic  algorithms  by  other  fire  control 
functions  are  not  violated.  A  numerical  integration  scheme  mi^t  be  made  to  compute 
time-of-fli^t  and  range  more  rapidly  tiian  the  China  Lake  Algorithm  since  the  weapon 
trajectory  is  relatively  strai^t  and  ranges  are  generally  small. 
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3. 2«  2  Hre  Control  System  For  Gunnery 


The  problem  of  formvilating  a  fire  control  solution  for  gunnery  modes  (and  also 
for  bombing)  involves  predicting  the  target  position  one  time-of-flig^t  (TOF)  in  the 
future  and  simultaneously  computing  the  weapon  ballistic  trajectory  out  to  one  TOF. 
For  consistency  the  estimation/prediction  problem  and  the  ballistic  trajectory  prob¬ 
lem  will  be  formulated  in  air  mass  coordinates. 

Various  vectors  ;^hich  enter  into  the  fire  control  solution  for  the  (AAG,  AGG) 
gunnery  modes  are  illustrated  in  Figure  11.  The  projectile  ballistic  vector  d 
is  broken  down  into  components  GA  =  and  aT  =  such  that  Rg  is  along  the 
total  initial  air  mass  velocity  vector,  ,  and  G  is  along  the  local  vertical, 

where  Vm  is  the  muzzle  velocity  vector  (projectile  velocity  relative  to  attacker)  and 
^  is  the  gun  velocity  relative  to  the  air  mass.  Both  vectors  and  ^  include  the 
effect  of  drag  which  tends  to  oppose  the  air  velocity  vector  along  the  projectile  tra¬ 
jectory.  While  the  projectile  travels  from  point  G  to  point  I  in  air  mass  coordinates, 
the  target  moves  from  point  T  to  point  F,  where  the  vector  TF  =  KL  is  the  kinematic 
lead  vector  relative  to  the  air  mass.  A  hit  occurs  when  point  I  coincides  with  point  F, 
or  when  the  vectors 'gT  and  GF  become  equal,  that  is; 

GI  =  Rg  +  G  =  -PQ  +  Pg  +  R  +  KL=^ 

This  is  the  fundamental  relationship  used  in  the  derivation  of  the  fire  control  solution 
and  the  aircraft  steering  command  signals,  based  on  Figure  11.  This  expression  is 
valid  for  both  tiie  AAG  and  AGG  modes  and  provides  commonality  for  the  gunnery 
modes.  Note  that  in  the  AGG  mode  the  target  actually  travels  relative  to  earth,  how¬ 
ever  for  consistency  KL  represents  the  motion  relative  to  the  air  mass. 
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Points: 


C 

S 

G 

T 

F 

A 

I 


Attacker  IMU  or  CG 
Attacker  Sight 
Attacker  Gun 
Present  target  position 

Future  target  peeUlon  1„  air  m.s  toordlnatea 
Aim  point 

Weapon  impact  point 


Vectors 

Pg  :  Parallax,  attacker  IMU  to  sight 

Pq  :  Parallax,  attacker  IMU  to  gun 

R  :  Present  range  to  target 

^  :  Kinematic  lead  vector  relative  to  the  air  mass 

G  :  Gravity  drop 

Rg  .  Weapon  travel  along  total  velocity  vector  relative  to  the  air  mass  V^+V^ 

Vq  :  Attacker  velocity  relative  to  the  air  mass  at  the  gun  station 

:  Muezle  velocity  along  gunline  (projectile  velocity  relative  to  attacker) 

FIGURE  11.  VECTORS  ASSOCIATED  WITH  AIR-TO-AIR  HRE  CONTROL  MODE 
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Future  Target  Position  relative  to  sight  (In  air  mass  coordinates) 


SF  -  R  +  KL 

Predicted  Weapon  Impact  Point  relative  to  sight  (In  air  mass  coordinates) 

SI  .Rg  +  G-f  Pg-Pg 

Predicted  Miss  Distance  Is  given  by  AM  relative  to  sight 


AM  = 


+ 


(3-5) 


(3-6) 


(3-7) 


where  AM, ,  AM_  and  AM_  are  the  components  of  vector  aM"  In  arbitrary  coordinates 

i.  A  O 


AM  = 


rAM^l 

AM, 


L  AM^J 


-P  +  P  +  R  +  ICL  -  G  -  R^ 


=  -  SI 


(3-8) 


Gunllne  Direction  Cosines  for  a  Hit 

The  orientation  of  the  gunllne  Is  defined  by  the  muzzle  velocity  vector 
Defining  the  distance  loss  due  to  drag  as  DB -  t^  (Vj^  +  Vq)  and  solving 
for  Vj^tj  from  Equation  3-4  yields: 


Gl  -  -  -Pj  +  Pj  +  E  +  Kl  -  DB,  -  G  - 

The  gunllne  direction  cosines  are  then  given  by 


(3-9) 


'gLj^/  I^I 
GL^/  1^1 
GL3/|^i 

where  IGLI  Is  the  norm  and  GL^^,  GL^,  and  GLg 
In  arbitrary  coordinates. 


(3-10) 


are  the  components  of  the  vector  GL 


Desired  Range  for  a  Hit  Is  by  definition  the  current  range  which  satisfies  the  condi¬ 
tion  for  a  hit  In  Equation  3-4.  Solving  for  R  and  designating  the  result  by  R^  yields: 

Rjj  =  R^  +  G  +  Pp  -  Pg  -  (3-11) 
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Gun  Lead  Angle  is  the  angle  between  the  desired  LOS  (or  vector 'Rjj)  and  the  gunline 
vector  Vj^tj.  The  direction  cosines  of  this  angle  can  be  obtained  by  resolving  the  gun- 
line  vector  Vj^tj  in  the  LOS  or  HUD  coordinates  (GL^,  GLg,  GL^)  and  computing  ihe 
direction  cosines  as  in  Equation  (3-10).  As  an  alternative  procedure,  the  azimuth 
(4T)  and  elevation  (dE)  angles  of  the  gun  lead  angle  relative  to  arbitrary  (body,  gun, 
LOS  or  HUD)  coordinates  can  be  obtained  as  follows: 


AT  =  tan 

AE  “  -sin 


+  sin 


(3-12) 


where  R-  ,  R^,  R-^  are  the  components  of  the  desired  range  vector  R, 
1  ^2  ^3 

GL^,  GLg  the  components  of  the  gunline  In  arbitrary  coordinates. 


and  GLj^, 


3. 2. 3  Tlme-of-Fllght  Computation  for  AAG  and  AGG 

The  purpose  of  the  fire  control  solution  for  AAG  and  AGG  is  to  obtain  the  weapon 
TOF  and  the  aiming  vector  direction  cosines  that  satisfy  the  condition  for  a  hit  In 
Equation  (3-4).  The  solution  method  for  the  TOF  computation  depends  on  whether 
closed-form  expressions  or  numerical  Integration  methods  are  used  for  the  weapon 
ballistics.  The  functional  requirements  for  each  case  will  be  obtained  below. 


3.2.3. 1  TOF  Computation  Using  Closed-Form  Ballistic  Equations 

The  TOF  Is  computed  from  the  scalar  equation  resulting  from  the  projection 
of  Equation  3-4  onto  a  convenient  axis,  such  as  the  LOS.  Since  the  vectors  Rg,  G 
and  KL  depend  on  the  TOF,  t^  must  be  solved  Iteratively  from  the  resulting  equations 
until  the  error  e  becomes  arbitrarily  small. 

'  -  S  -  "'Sj  -  '*1  -  '‘'-1  (3-13) 

where  subscript  (1)  denotes  the  vector  component  along  the  arbitrary  axis.  The 
Newton-Raphson  algorithm  for  solving  t^  Is; 
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where 


(3-14) 


3Gj^  3KLj^ 

atj  ■'■  3tj  “  3t^‘ 


3 . 2. 3. 2  TOF  Compulation  Using  Numerical  Integration  Methods 

China  Lake  algorithms  can  be  modified  to  converge  on  slant  range  Instead  of 
altitude  for  the  numerical  Integration  of  weapon  ballistic  equations.  This  procedure 
would  be  equivalent  to  the  closed-form  Iteration  along  the  LOS  proposed  In  Section 

3.2.3. 1. 


The  modified  China  Lake  algorithm  would  compute  the  weapon  TOF  and  the 
weapon  ground  range  from  given  Initial  velocity  and  slant  range,  as  shown 

In  Figfure  12.  From  the  geometry  in  this  figure  the  vectors  Rg  and  G  used  in  the 
fire  control  solution  are  obtained  as  follows: 


cos 


2  2  — 
Rg  -  Rj,  +  Rg  tany)  Z 


(3-15) 


G  *  (/ 

where  Z  Is  a  imit  vector  along  the  local  vertical,  Rg  is  slant  range  and  y  is  the  eleva¬ 
tion  angle  of  the  total  velocity  vector,  +  Vq  measured  positive  counterclockwise. 


FIGURE  12.  MODIFIED  CHINA  LAKE  NUMERICAL 
INTEGRATION  ALGORITHM  FOR  AAG  AND  AGG 
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An  iterative  solution  can  be  set  up  to  adapt  the  China  Lake  algorithm  to  FIREFLY 
gunnery  solutions  as  shown  In  Figure  13.  In  this  diagram  initial  values  for  slant 
range  (Rg)  and  velocity  vector  (Vj^  +  V^)  are  assumed.  Weapon  ballistic  equations 
are  integrated  numerically  to  obtain  t,^  and  Rq.  Knowing  t^  the  kinematic  lead  vector 
and  the  miss  distance  are  computed,  and  the  Initial  slant  range  Rg  is  adjusted  acc.ord- 
ingly  In  an  iterative  manner  imtll  the  miss  distance  Is  reduced  to  an  acceptable  level. 


Rg.  G 


FIGURE  13.  time  OF  FLIGHT  COMPUTATION  USING  A 
MODIFIED  (CHINA  LAKE)  NUMERICAL  INTEGRATION  ALGORITHM 


28 


3.2.4  Fire  Control  System  for  Bombing 

FIREFLY  bombing  algorithms  albw  for  a  maneuvering  approach  to  bomb  re¬ 
lease  in  order  to  enhance  survivability  to  enemy  (AAA  or  SAM)  defenses  while 
maintaining  a  high  level  of  offensive  weapon  delivery  accuracy.  The  maneuver  con¬ 
sists  of  a  circular  flight  trajectory  on  a  nonhorizontal  turning  plane.  The  release 
point  along  this  trajectory  is  computed  continuously  (CCRP),  by  assuming  constant 
turn  rate  and  constant  ownship  velocity  relative  to  the  air  mass. 

Figure  14  shows  the  turning  plane  geometry  for  bombing.  ’  The  turn  is 
initiated  at  point  C  and  the  bomb  is  released  at  point  R.  Point  C  is  also  the  point 
where  the  target  is  sighted,  and  the  CCRP  computation  is  performed.  From  point  C 
to  point  R  the  aircraft  flies  a  circular  trajectory  relative  to  the  air  mass.  The  veloc¬ 
ity  vectors  relative  to  the  air  mass  at  points  C  and  R  are  tangent  to  the  circle;  hence, 
the  following  fundamental  relationship  holds  for  FIREFLY  bomblno’. 

2  2 

„  (3-16) 

CA 

In  vector  form  this  relationship  beco’ues 


2  2 

S  2V_.  “ 


(3-17) 


The  vector  diagram  used  for  computing  the  desired  weapon  Impact  point  (F)  is  shown 
in  Figure  15.  From  this  figure  the  vector  CF  is  given  by: 

S  of 


The  vector  diagram  at  the  release  point  is  shown  in  Figure  16,  wherein  the 
combined  effect  of  bomb  ejection  velocity  (V^,),  wind  velocity  W  and  parallax  (Pg  + 
w  X  Pg  t^)  is  represented  as  an  impact  point  correction  (PI)  in  order  to  allow  the  funda¬ 
mental  relationship  in  Equations  (3-16  and  17)  to  be  used  for  the  generalized  bombing 
solution.  Otherwise  a  point  P,  directly  above  the  weapon  impact  point  on  the  turning 
plane  which  lies  along  the  velocity  vector  cannot  be  defined.  In  general,  however, 
the  incremental  bomb  initial  velocity  Vg  +  w  x  P^  will  affect  the  aerodynamic  forces 
acting  on  the  bomb.  But,  this  effect  is  negligibly  small. 
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CIRCULAR  FLIGHT  PATH 


R 

P 


IMU  location  (or  aircraft  CG) 
at  turn  initiation 

Bomb  release  point 

Point  on  turning  plane 
directly  above  weapon 
impact  point 

Weapon  impact  point 


"V 

G’ 

Rn 


Unit  vector  along 

Vector  from  point  C  to  point  P 

Weapon  travel  along  from 
release  to  point  ?. 

Vertical  drop  of  weapon 

Radius  of  turning  circle 


X  Angle  between  V^^  and  Rp 
FIGURE  14.  TURNING  PLANE  GEOMETRY  FOR  BOMBING 


POINTS 
C 
S 
T 
F' 


IMU  location  at  turn  initiation 

Target  tracking  sensor 

Present  target  position 

Future  target  position  relative 
to  the  air  mass 

Desired  weapon  impact  point 
relative  to  the  air  mass 


KL  Kinematic  lead  vector  relative 

—  to  the  air  mass 

of  Target  position  offset 


FIGURE  15.  VECTOR  DIAGRAM  AT  TURN  INITIATION 
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VECTORS 


W  Wind  velocity  relative  to 
inertial  coordinates 

?_  Bomb  parallax 

D 


Vg  Ownship  velocity  at  bomb 
location 


Vg  Bomb  ejection  velocity 
relative  to  body 


Defined  above. 


FIGURE  16.  VECTOR  DIAGRAM  FOR  BOMB  RELEASE  POINT  COMPUTATION 


The  condition  for  a  hit  is  satisfied  when  the  desired  impact  point  (point  F  in 
Figure  15  coincides  with  the  actual  impact  point  (point  I  in  Figure  16.)  From  this 
condition  the  vector  CP  =  Rp  defined  in  Figure  3-12  is  obtained  as  follows. 


Rp  =•  Pg  +  R  +  KL  +  -  W  tg  -  (Vg  +  u)  X  Pg)  '  t^  -  G'  -  Pg*  (3-19) 

The  angle  subtended  by  the  circular  path  CR  is  wt^,  hence  from  the  geometry  in 
Figure  14  the  time  to  go  before  release,  t  ,  becomes: 

O 


t 

g 


(3-20) 


where  |Rp  -  Rp|  denotes  the  length  of  vector  CR. 


Fire  Control  Parameters 

Various  parameters  which  are  needed  for  the  bombing  solution  can  now  be 
expressed,  as  a  .function  of  the  vectors  defined  above. 


(*)  Note  that  if  KL  is  relative  to  inertial  coordinates  Wt,  should  be  substituted  in  place 


of  Wt 


g' 
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Weapon  Release  Point,  or  distance  to  go  to  bomb  release 


(3-21) 


Future  (offset)  Target  Position  relative  to  tracking  sensor  location  S, 

^=R+^+D, 

of 


(3-22) 


Predicted  Weapon  Impact  Point  relative  to  point  S. 

sl  =  -  Pg  +  G'  +  (Vg  +  u  X  Pg)  •  +  Pg  +  wt^ 


(3-23) 


Predicted  Mias  Distance  Is 


m 


= 


2  2  2 
+  AM2  +  AM^ 


(3-24) 


where  /IM2,  -dMg  are  the  components  of  the  vector  AM  in  arbitrary  coordinates, 


AM  =•  SF  -  SI  »  R  +  KL  +  -  Rp  +  Pg  -  G'  -  Wtg 

-  (V^  +  :  X  P3)  •  -  p. 


(3-25) 


Condition  for  a  Hit  is 


AM  =  0  or  SF  =  SI 


(3-26) 


Time  to  Bomb  Release  Command 

The  equations  above  are  derived  relative  to  the  true  release  point  R.  Because 
of  the  rack  delay  (t)  the  bomb  release  must  be  initialed  r  seconds  before  reaching 
point  R  or  t  -  r  seconds  after  passing  through  point  C,  hence 

o 


Time  to  go  2  _i 

before  bomb  release  =  —  sin 
command 


1^  ~  ^1 
2  V 


-  T 


(3-27) 


3.  2. 5  Tlme-of- Flight  Computation  for  the  Bombing  Solution 

The  vectors  defined  above  depend  on  the  TOF  and  hence  the  TOF  must  be  com¬ 
puted  in  order  to  complete  the  fire  control  solutions  and  determine  the  release  point 
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along  the  turning  circle.  :..ie  TOF  and  Qie  control  laws  are  derived  from  the  funda¬ 
mental  relationship  for  FIREFLY  bombing  in  Equation  3-17.  The  control  law  adjusts 
the  inclination  of  the  turning  plane  (the  direction  of  u)  v/hile  the  TOF  is  computed  such 
that  the  magnitudes  of  the  vectors  in  Equation  S-17  are  made  equal. 

3. 2. 5. 1  TOF  Computation  Using  Closed- Form  Balliatics 

The  magnitude  equality  In  Equation  3-17  is  satisfied  b}^  nulling  the  error  e 

‘  ""  ®D2  "  ®D1 

where 

^2-^Dl 

*D1  ■  ^  [V  \  -  V  *  (3-28) 

I<B2-  |S,xi4,i  -  |Cj| 

Vector  notation  is  used  to  derive  expressions  for  the  partial  derivatives  of  R^  and  Rj^ 
with  respect  to  t^.  Hence,  the  formulation  presented  here  is  more  general  than  the 
Firefly  II  formulation  and  does  not  require  simplifying  assumptions,  similar  to  those 
in  the  FIREFLY  II  report,  to  be  made  in  order  to  obtain  the  desired  expressions. 

All  quantities  in  Equation  3-28  are  nonlinear  functions  of  the  TOF,  hence  must 
be  solved  iteratively  from  these  equations.  The  Newton-Raphson  algorithm  for  solving 
t^ is  given 


where 


(3-29) 


(*)  Superscript  T  denotes  vector  transposition. 
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3.  2. 5.  2  TOF  Computation  Using  Numerical  Integration 

A  particalarlj'  useful  method  for  computing  the  weapon  ballistic  trajectory  is  to 
integrate  the  ballistic  equations  numerically  from  bomb  release  to  impact  by  using 
the  China  Lake  algorithms.  These  algorithms  provide  better  weapon  delivery  accuracy, 
but  are  not  directly  amenable  to  the  FIREFLY  bomb  algorithm.  The  China  Lake  algo¬ 
rithms  compute  the  TOF  and  ground  range  (Rq)  by  specifying  the  initial  bomb  release 
altitude  (h)  (above  impact  point)  and  the  initial  vertical  and  horizontal 

components  of  total  weapon  velocity  at  bomb  release,  as  illustrated  In  Figure  17. 


=  Ground  range 


■=  TOF 


Weapon  altitude 


Velocity 
relative  to  the 
air  mass 


FIGURE  17.  CHINA  LAKE  NUMERICAL  INTEGRATION  ALGORITHMS 


From  the  geometry  in  Figure  17  and  from  the  definitions  given  above  the 
vectors  Rj^  and  G*  which  define  the  bomb  ballistic  trajectory  are  obtained  as  follows; 


G' 


CA 


CA 


cosy 


(3-30) 


(h  +  Rq  tan  y)  Z 


where  Z  is  unit  vector  along  the  local  vertical  and  y  is  the  elevation  angle  of  VCA 
relative  to  the  horizontal  plane  (measured  positive  counterclockwise). 

An  iterative  solution  similar  to  the  one  used  for  the  AAG  and  AGG  must  be 
developed  if  numerical  integration  is  to  be  used  effectively  in  FIREFLY  bombing. 
This  solution  can  be  made  to  converge  on  altitude  or  on  slant  range  depending  on 
which  method  gives  better  results. 
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3, 2. 6  Fire  Control  Solution  for  Stick  Bombing 


The  turning  rate  for  FIREFLY/slngle  target  bombing  is  an  independent  variable 
which  can  be  chosen  conveniently  by  the  pilot.  This  is,  the  pilot  has  the  option  of  ini¬ 
tiating  the  turn  immediately  or  continue  at  the  present  course  and  then  make  a  tighter 
turn.  For  stick  bombing  the  turning  rate  is  no  longer  an  independent  variable  at  the 
pilot's  disposal;  it  must  be  computed  along  with  the  point  at  which  the  turn  is  initiated 
by  the  FIREFLY/Stlck  bombing  algorithm. 

The  geometry  for  horizontal  stick  bombing  is  shown  in  Figure  18.  Since  the 
altitudes  and  ground  ranges  associated  with  the  end  points  of  the  target  are  equal,  there 
exists  a  unique  release  circle  with  Center  C  along  the  perpendicular  bisector  AT  for 
which  its  radius  satisfies  the  relationship  in  Equation  3-16.  The  release  circle  is 
tangent  to  the  velocity  vector  at  point  B,  and  hence,  the  flight  trajectory  consists  of  a 
straight  line  segment  AB  followed  by  the  circular  path  BR^  R^.  The  stick  bombing 
algorithm  will  be  extended  to  the  nonhorizontal  turning  plane  case  in  the  GM  design 
phase. 

The  weapons  will  be  released  between  points  R-  and  Rg  on  the  turning  circle 
where  the  segments  R^^  and  Rg  are  tangent  to  the  circle.  The  distance  from  any 
point  along  the  linear  flight  segment  AB  to  the  turn  initiation  point  B  is  referred  to  as 
"the  distance  to  go. "  In  addition  to  this  distance,  the  fire  control  solution  for  stick 
bombing  must  also  compute 

•  The  turning  rate,  and 

•  The  release  points  along  the  circular  trajectory. 


URE  18.  FIREFLY /STICK  BOMBING  GEOMETRY  FOR  HORIZONTAL  TURNING  PLANE 


3. 3  STATE  ESTIMATION  AND  FUTURE  TARGET  POSITION  PREDICTION 


State  estimators  compute  the  required  states  of  owp.ship,  atmosphere,  and  target . 
using  subsystem  inputs  and  reasonable  dynamic  equations.  The  target  state  estimator 
shall  be  common  for  both  aerial  and  ground  targets;  only  the  data  constraints  will  be 
different.  The  future  target  position  and  the  kinematic  lead  vector  shall  be  computed 
using  target  state  estimates  and  the  weapon  TOF. 

3.3.1  Ownship  State  Estimator  Requirements 

Implementation  of  the  FIREFLY  fire  control  solution  requires  knowledge  of  the 
ownship  states  with  which  to  reference  the  target  relative  motion  and  obtain  the  we^on 
initial  conditions.  Some  ownship  states  such  as  angular  rates  and  aircraft  orientation 
are  measured  directly  (or  processed)  by  avionics  sensors.  To  differentiate  from  these 
states,  only  the  ownship  states  which  are  generated  through  an  estimation,  filtering,  or 
data  mixing  function  are  considered  a  part  of  the  ownship  state  estimation  process. 
These  include; 

•  Ownship  accelerations  (A_) 

•  Ownship  inertial  velocity  (V^) 

•  Angle-of-attack  (a)  and  sideslip  (fi) 

3. 3. 1.1  Velocity  and  Acceleration  Estimators 

Inertial  axes  are  considered  to  be  fixed  rotationally  in  inertial  space  with  the 
origin  at  the  center  of  the  earth  for  computations  involving  inertial  quantities. 

Inertial  computations  incorporating  doppler  measurements  shall  consider  the 
effects  of  earth  curvature  as  a  fimction  of  altitude  if  applicable.  Ownship  velocity  and 
acceleration  are  computed  by  one  or  more  of  the  following  sensor  outputs. 

Stabilized  Platform;  Sensed  accelerations  from  the  INS,  (Ag)  can  be  combined 
with  the  gravitational  acceleration  (A^)  to  provide  an  estimate  of  total  body  accelera¬ 
tion  resolved  in  appropriate  coordinates  (such  as  body)  as  shown  below. 


[bp]  -  platform  to  A/C  body  transformation 
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Strapdown  -  Sensed  accelerations  in  body  coordinates  can  be  summed  with  the  gravity 
vector  resolved  through  the  HARS,  Vertical  gyro  attitude  angles  to  estimate  the  total 
ownship  acceleration  as  in  Equation  3-32. 


*C2-*S*  ™  *0 


(3-32) 


where  [BH]  is  HARS  axes  to  body  axes  transformation. 

Similarly,  ownship  inertial  velocity  can  be  obtained  directly  from  the  INS  measure¬ 
ment.  If  such  measurements  are  not  available,  ownship  inertial  velocity  can  be  calcu¬ 
lated  from  doppler  radai’  velocity  measurements  (or  estimates)  Vd  as  follows: 


''c  - ’d  +  “  * 


(3-33) 


where  Q  is  Earth’s  angular  rate,  R^,  is  Earth's  radius  and  x  R^  is  calculated  in  the 
fire  control  computer. 

An  alternate  means  cf  computing  inertial  ownship  velocity  involves  fusing  and  in¬ 
tegrating  either  or  both  stabilized  platform  and  strapdown  accelerometer  outputs, 

A„  ,  A-  according  to 
^2 


V^(i)  =  V„(i-i)  +r  Ik  A„  (T)  +  (1-K)  A  (T)"]  dT  (3-34) 

^  ^  ^  (i-i)A  '^2  J 


where  K  is  a  coefficient  ranging  from  0  to  unity  and  d  is  the  update  time. 
g,  B  Estimator 

The  direction  cosines  (a^,  of  the  air  velocity  vector  relative  to  ownship 
body  coordinates  and  the  associated  true  Euler  angles  (  q  and  p )  are  given  by 


01*  =  V  /V 

caw'  '^tas 

3^ 

=  V  /V 

CAV'  TAS 

a  =  tan“^  /^CAW  ]  ; 

3 

•=  Sin"^  ('^Ckv\ 

'^"CAU^ 

Vvtas/ 

(3-35) 


*Vtas  islhe  air  speed  at  the  IMU  or  aircraft  CG,  obtained  by  applying  the  proper 
corrections  to  measured  air  speed. 
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These  equations  can  be  used  to  compute  the  direction  cosines  and  the  true  Euler  angles 

directly  from  the  air  velocity  vector  An  alternate  method  for  estimating  the 

1  1  1  1 
direction  cosines  a  and  is  shown  in  Figure  19.  In  this  figure  a  and  4  are 

estimated  first  from  acceleration  (A^,  A^)  and  vane  (o^,  measurements  through 

a  Kalman  filter  with  constant  gains  Kg,  Kg,  and  K^.  Integrating  the  and 

estimates  then  yields  the  desired  estimates  of  the  direction  cosines  and  4^. 

3. 3.  2  Atmosphere  Estimators 

The  specific  atmospheric  parameters  that  must  be  measured  or  estimated  in¬ 
clude  ownship  velocity  relative  to  the  air  mass,  air  pressure,  wind  velocity  (both 
horizontal  and  vertical),  air  density,  and  air  pressure. 

Table  2  is  a  summary  of  the  assumptions  made  concerning  some  of  the  atmos¬ 
pheric  parameters.  These  assumptions  were  made  for  the  atmospheric  parameter/ 
weapon  delivery  mode  denoted  by  "No  Requirement. "  hi  addition,  the  table  also  denotes 
for  which  parameters  there  is  a  requirement  to  be  measured  or  estimated  for  each 
weapon  delivery  mode.  The  numbered  keys  in  the  boxes  refer  to  the  reasons  corre¬ 
sponding  to  the  assumptions  made  or  requirements  established,  and  are  further  dis¬ 
cussed  in  the  following  section. 


TABLE  2 

ATMOSPHERIC  PARAMETER  REQUIREMENTS 


Weapon 

Delivery 

Air 

Velocity 

Wind 

Air 

Mode 

Horizontal 

Vertical 

Shear 

Density 

a/a  Guns 

Yes 

No  Requirement 

No  Requirement 

No  Requirement 

Yes 

(1) 

(1) 

(1).  (3) 

(4) 

A/G  Guns 

Yes 

Yes 

(2) 

Yes 

(2) 

No  Requirement 
(3) 

Yes 

(4) 

Bombing 

Yes 

Yes 

(2) 

Yes 

(2) 

No  Requirement 
(3) 

Yes 

(4) 

Reasons  for  assumptions/requirements: 

(1)  All  wind  is  ignored  for  AAG. 

(2)  Bias  due  to  uncompensated  wind  is  too  large  to  ignore. 

(3)  No  way  to  measure  wind  shear. 

(4)  Required  for  all  air  data  measurements. 
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For  the  AAG  mode  it  is  reasonable  to  ignore  all  wind  and  wind  shear  effects. 

If  the  air  mass  at  the  ownship  and  target  positions  (and  the  space  in  between)  is  moving 
at  a  uniform  rate  with  respect  to  the  ground  (i.  e. ,  no  wind  shear),  then  for  all  practical 
purposes  the  wind  velocity  can  take  on  any  value  at  all  and  it  would  not  matter.  What 
matters  is  target  motion  relative  to  ownship. 

For  the  AGG  mode,  the  average  air  mass  velocity  relative  to  the  ground  will 
affect  the  fire  control  solution  since  the  aircraft  moves  relative  to  the  air  mass  while 
the  target  moves  relative  to  the  ground.  However,  if  the  bullet  trajectory  computation 
and  future  target  position  prediction  are  both  performed  in  air  mass  coordinates,  wind 
velocity  will  not  enter  into  the  computation. 

In  FIREFLY  bombing,  the  taiget  acquisition  and  the  CCRP  computation  point 
differs  from  the  release  point.  Since  the  aircraft  moves  relative  to  the  air  mass,  the 
actual  release  point  would  move  by  Wt^  (W  =  wind  velocity,  t^  =  time  to  go  from  ac¬ 
quisition  to  bomb  release)  relative  to  the  release  point  predicted  in  the  air  mass. 
Therefore,  wind  velocity  must  be  estimated  for  tlie  bombing  niode.  Wind  will  also 
affect  the  bomb  trajectory  througli  air  speed  which  produces  drag. 

At  first  appearance  it  might  seer  ‘^asonable  to  ignore  all  vertical  wind  since 
vertical  motion  of  the  atmosphere  is  one  to  three  orders  of  magnitude  smaller  than 
horizontal  wind.  However,  even  though  the  actual  vertical  wind  velocity  is  small,  the 
contribution  of  neglecting  the  vertical  wind  to  the  total  weapon  miss  distance  is  not 
necessarily  small.  To  illustrate  this.  Table  3  gives  the  ballistic  range  sensitivity 
to  an  error  in  the  vertical  and  horizontal  wind  velocities  for  five  different  release  con¬ 
ditions  and  for  a  low  drag  bomb.  These  sensitivities  were  generated  using  Northrop's 
air-to-ground  weapon  delivery  trajectory  simulation  program. 

TABLE  3.  SENSTWiry  OF  WEAPON  IIvIPACT  POINT  TO  VERTICAL 

AND  HORIZONTAL  WIND 


Release  Condition 


Dive  Angle  (deg.) 
Altitude  (ft.) 


Vertical  Wind 

Sensltivlty(ft/ft/sec)  21.648 


Horizontal  Wind 
Sensitlvity(ft/ft/sec)  7.383 


11.117  6.745 


6.379  6.112 


60 

10,000 


5.001  4.102 


7.463  11.645 
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The  air  velocity  vector  is  computed  from  true  air  speed  measurement 
Vtas  heading  information.  Initially,  wind  direction  and  magnitude  information 
is  received  from  the  meteorology  station.  Knowing  the  wind  vector  and  the  ground 
speed  from  the  INS  gives  the  air  velocity  vector  VCA  which  closes  the  wind  triangle. 
Subsequently,  a  corrected  wind  velocity  is  computed  by  simply  subtracting  air  velocity 

A 

vector  estimate  ^^0“^  ownshlp  inertial  velocity  estimate  V^. 


W  - 


(3-36) 


3.3.3  Target  State  Estimator  Requirements 

The  primary  function  of  the  target  state  estimator  is  to  provide  accurate  estimates 
of  target  position,  velocity,  and  acceleration  by  processing  ownship  states  and  desig¬ 
nated  target  position  Information.  Target  position  Information  is  either  supplied  by 
the  pilot  as  a  position  offset  in  inertial  coordinates  or  measured  through  tracking 
sensors  such  as  EO  sensor,  radar,  and/or  laser. 

The  important  ingredients  which  affect  the  target  state  estimator  structure  are 
the  coordinate  frames  chosen  for  the  estimation  process,  the  target  model,  and  the 
observations. 

3.3. 3.1  Coordinate  Frames 

Possible  coordinates  for  the  estimation  process  include  inertial,  ownshlp  lx)dy, 
and  polar  LOS  coordinates.  The  (slm)  LOS  coordinates  are  defined  in  Section  4. 5 
where  &  is  along  the  LOS  and  I  and  m  are  perpendicular  to  it.  Ownshlp  body  coordinates 
can  be  eliminated  based  on  the  argument  that  the  accuracy  of  the  estimates  will  be 
affected  adversely  by  body  rotation  as  compared  with  inertial  and  polar  LOS  coordinates. 

Both  inertial  and  polar  LOS  coordinates  have  been  used  successfully  in  target 
state  estimation;  the  F-16  uses  inertial  coordinates  whereas  the  F-15  and  the  FIREFLY 
II  report  use  of  polar  LOS  coordinates.  There  are  distinct  advantages  associated  with 
each  of  these  coordinate  frames. 

Target  motion  can  be  best  modeled  using  inertial  coordinates  (fixed  relative  to 
earth  or  relative  to  the  air  mass),  while  the  LOS  measurements  can  be  best  defined 
using  the  polar- LOS  coordinates.  Therefore,  the  choice  of  an  appropi’iate  coordinate 
system  for  the  generalized  target  state  estimator  design  should  be  based  on  a  trade-off 
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study  to  be  conducted  as  part  of  the  design  effort.  Usually,  the  choice  of  coordinate 
frames  is  dependent  upon  the  sensor  characteristics. 

Generalized  target  state  vectors  compatible  with  inertial  and  polar  LOS  coordi¬ 
nates  are  given  by; 


(Inertial) 


\ ' 

R 

R 

R 

y 

R 

z 

"s 

^RX 

(Polar-LOS) 

^2  “ 

T 

S 

V 

RZ 

tm 

a^y 

(3-37) 


The  components  of  the  state  vector  Xj  include  target  range  R,  relative  Vv..o  y  Vj^, 
and  target  acceleration  resolved  in  (XYZ)  inertial  coordinates.  The  components  of 
Xg  include  range  (R),  range  rate  (R),  elevation  (Eg)>  and  azimuth  (T^)  angles  of  the 
LOS  measured  relative  to  ownship  body  coordinates,  angular  rates  of  the  LOS  perpen¬ 
dicular  to  the  LOS  (w. U4._),  and  the  target  acceleration  resolved  in  rotating  (siro.) 
sight  coordinates. 

Estimates  in  either  inertial  or  polar  LOS  coordinates  are  transformed  to  ownship 
body  coordinates  by  using  the  appropriate  transformation  naatrices.  The  (s^m)  coordi¬ 
nates  are  roll-stabilized  in  which  no  rolling  motion  takes  place  about  the  s-axis. 

To  retain  flexibility,  the  generalized  state  estimator  could  be  designed  to  operate 
with  both  sets  of  cooi'dinates  (inertial  and  polar  -LOS)  by  closing  the  appropriate 
switches  shown  in  Figure  20).  The  state  transition  and  propagation  matrices  0,  G, 
and  r  are  invariant  with  the  set  of  sensors  used  on  a  specific  aircraft  and  depend  only 
on  the  coordinate  frame  selection.  The  observation  matrix  H,  on  the  other  hand, 
depends  on  both  tlie  coordinate  frame  and  measurement  complement.  In  this  way  the 
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FIGURIS  20.  GENERALIZED  TARGET  STATE  ESTIMATOR 


impact  of  sensor  availability  on  estimator  design  is  localized  to  the  selection  of  an 
appropriate  matrix  H,  and  the  estimator  is  thus  modularized.  In  addition  to  the  ob¬ 
servation  vector  Y  and  the  measurement  error  covariance  matrix  L  ,  the  ownship  ac¬ 
celeration  at  the  tracking  sensor  Ag  must  also  be  supplied  as  inputs  to  the  target  state 
estimator.  The  latter  is  obtained  from  ownship  acceleration  at  the  IMU,  A^,  and  body 
rates  was  indicated  in  Figure  20.  Target  state  estimates  are  converted  from  LOS 
or  inertial  coordinates  to  body  coordinates  where  the  fire  control  solution  is  computed. 


3.3. 3. 2  Target  Model 

A  representative  dynamic  model  which  defines  the  position  (P.p),  velocity  (V.p), 
and  acceleration  (A.p)  of  the  target  relative  to  inertial  coordinates  at  any  given  time  is: 


d 

dt 


0  10 

0 

= 

0  0  1 

+ 

0 

A 

0  -Y  -e 

A 

u 

(3-38) 


where  a  random  (white)  process,  u,  is  added  to  the  model  to  introduce  randomness 
into  the  target  motion.  0  and  y  are  appropriate  3x3  matrices  which  define  time  cor¬ 
relation  and  interaction  effects.  This  target  acceleration  model  is  second-order 
Gauss-Markov  and  has  sufficient  flexibility  to  allow  a  wide  variety  of  targets  with 
varying  degrees  of  time  correlation  to  be  represented  by  proper  choice  of  the  0  and  y 
matrices.  For  instance,  setting 0=0  will  produce  a  circular  target  trajectory  for 
modeling  an  evasive  turn  in  the  AAG  mode.  Similarly,  setting  y  =  0  would  allow 
various  ground  targets  to  be  represented;  0=0  corresponds  to  constant  acceleration 
and  ;9=  00  to  a  rapidly  changing,  unpredictable  acceleration  pattern. 

3. 3. 3. 3  Observations 

Observations  will  usually  involve  range  to  target,  LOS  direction  cosines,  and 
in  most  cases  it  will  include  the  angular  rates  of  the  LOS.  These  observations  will 
be  obtained  from  the  measurement  function  discussed  in  Section  3.  5. 
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3.3.4  Kinematic  Lead  Vector  and  Future  Target  Position 

The  target  dynamic  model  in  Equation  (3-38)  can  be  used  to  predict  the  target 
position  on  TOP  in  the  future.  In  general,  the  future  target  position  vector  is  given 
by 

P(tf)  -  •  P(0)  +  03^2^*=^  *  **13^*"^ 

where  ^3^3^ (tj),  ^3^2^*'f^'  ^13^*'f^  3x3  state  transition  matrices  related  to  the  model 

in  (3-38)  and  P(0),  V.j,(0),  and  A^(0)  are  current  estimates  on  target  position,  velocity, 

and  acceleration  obtained  from  the  target  state  estimator.  V„,  is  related  to  the  esti- 

A  T 

mated  relative  target  velocity  vector  Vj^  by 

*  'CA  +  “ 

The  kinematic  lead  vector,  KL,  is  defined  as  the  distance  traveled  by  the  target 
during  one  TOP  and  is  given  by 


KL  -  P(t^)  -  P(0) 


(3-41) 


A  familiar  expression  for  the  kinematic  lead  vector  can  be  obtained  by  setting 
P-y  ^  0  in  the  model  in  Equation  (3-41).  The  result  is; 


KL 


(3-42) 
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3.4  CONTROL  LAW 


The  primary  purpose  of  aircraft  control  laws  is  to  compute  the  desired  kinematic  • 
responses  of  the  aircraft  based  upon  the  fire  control  solutions  and  to  supply  controlled 
inputs  to  the  flight  control  system  (FCS). 

3.4.1  Control  Laws  for  AAG  and  AGO 

For  AAG  and  AGG  weapon  delivery  modes  the  aircraft  angular  rate  command 
inputs  are  derived  from  the  fundamental  FIREFLY  control  philosophy  which  is  "to  make 
the  aircraft  angular  velocity  vector  U  equal  to  the  predicted  angular  rate  associated 
with  the  hiture  target  position  and  simultaneously  null  the  angular  error  between  the 
present  range  vector  R  and  the  desired  range  vector  for  a  hit  R^. "  This  control 
philosophy  is  illustrated  in  Figure  21  below . 


(a)  First  Firefly  Requirement  for  a  Hit  (u^  =  w) 


Target 
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According  to  the  FIKEFLY  n  control  philoaophy,  pitch  and  roll  rate  commands 
provide  coarse  adjustment  in  the  aircraft  flight  path  to  meet  the  requirements  for  a  hit, 
while  yaw  rate  command  provides  fine  tuning  in  gun  aiming  to  achieve  a  hit.  The  angu¬ 
lar  error  between  the  desired  range  and  present  range  R  is  nulled  by  turning  the 
aircraft  until  Rjj  becomes  collnear  with  R. 


The  control  law  is  given  by 

‘Ic  -  ‘ICG  ^  « 
"c  -  "CG  \  « 


w 


(3-43) 


Ki  (q  •  •  r)  +  (q  •  -  r  •  e^) 


where  K^,  are  appropriate  feedback  gains,  M  is  a  variable  gain  function 

and  e  ,  e  ai 
V  w 

in  Figfure  22. 


and  e  ,  e  are  the  elevation  and  azimuth  error  angles.  The  control  system  is  shown 

Y  Vr 


To  achieve  the  desired  aircraft  response  to  angular  rate  commands  the  flight 
control  system  must  have  a  structure  similar  to  that  shown  in  Figure  23.  In  this 
figure  the  proportional  plus  integral  controller  in  the  forward  loop  insures  unity  trans¬ 
fer  function,  in  the  steady  state,  between  commanded  rates  and  achieved  rates. 


3.4.2  Additional  Control  Laws  for  the  AGG  Mode 

In  the  AGG  mode,  the  interaction  between  the  longitudinal  and  lateral  aircraft 
dynamics  must  be  minimized  to  achieve  the  desired  response  fidelity  to  angular  rate 
commands  while  performing  large  evasive  maneuvers.  Therefore,  the  control  laws 
for  AGG  should  provide  essentially  decoupled  aircraft  responses  to  pitch  rate  and  yaw 
rate  commands.  Decoupled  response  is  also  desirable  for  the  AAG  mode  but  not  neces¬ 
sary,  because  for  the  AAG  maneuvers  the  coupling  between  the  lateral  and  longitudinal 
airci'aft  d5mamice  is  small;  hence,  omission  of  decoupling  control  laws  will  not  cause 
severe  degradation  in  weapon  delivery  accuracy. 

Additional  feedback  control  laws  are  required  to  decouple  the  lateral  and  longi¬ 
tudinal  aircraft  responses  for  the  AGG  mode.  The  decoupling  control  lav/s  xised  in  the 
GE  FIREFLY  II  report  consist  of  feedbacks  to  cancel  out  the  inertial  and  the  aerody¬ 
namic  coupling  terms  in  the  aircraft  equaticais  of  motion.  These  control  laws  essentially 
provide  neutral  stability. 
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FIGURE  22.  AIRCRAFT  CONTROL  SYSTEM  FOR  AIR-TO-AIR  GUNNERY 


The  pitching  moment  and  yawing  moment  equations  of  the  decoupled  aircraft  have 
the  following  form  when  linearized  about  the  steady  flight  condition: 

q  -  .q 

(3-44) 


r  ”  p  +  N  p  +  N  r 
Izz  r 


where  Mq,  Np,  Nr  are  stability  derivatives  and  Ixz,  Izz  are  aircraft  inertias. 


Decoupling  control  laws  similar  to  fliose  given  in  the  GE  FIREFLY  II  report 
should  be  Implemented  in  the  flight  control  system  to  achieve  the  desired  decoupled 
aircraft  response. 


3.4.3  CONTROL  LAWS  FOR  BOMBING 


The  command  signals  that  steer  the  aircraft  automatically  to  the  release  point 
are  derived  from  the  fundamental  bombing  Equation  (3-17).  The  objective  of  the  con¬ 
trol  law  is  to  adjust  the  orientation  of  the  angular  rate  vector  by  controlling  aircraft 
roll  rate  until  it  becomes  colinear  with  the  vector  Cp. 


= 


(3-45) 


As  shown  in  the  GE  FIREFLY  II  report,  the  roll  rate  command  is  given  by 


where  is  the  component  of  the  eiror  vector  7 


(3-46) 


r  -  (%  K  »  (3-«) 

CA. 

in  u-body  axis  and  Kp  an  appropriately  chosen  feedback  gain. 

In  addition  to  generating  the  steering  commands  for  conventional  aircraft  the 
generalized  mechanization  for  bombing  must  also 

•  Provide  steering  commands  for  the  A- 3  0,  and 

•  Be  compatible  f^r  stick  bombing. 

In  the  FIREFLY  approach  the  pilot  sets  the  magnitude  of  cJ  by  controlling  turn 
rate,  the  direction  of  is  made  equal  to  the  direction  of  Cp  by  controlling  roll  rate 
and  the  TOF  is  computed  from  the  magnitude  equality  of  the  vectors  on  both  sides  of 
Equation  (3-17).  The  TOF  determines  the  release  point  relative  to  target  position. 
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Two  alternate  approaches  for  using  the  pitch  control  axis  of  the  A-10,  instead 
of  the  roll  control  axis,  for  steering  the  aircraft  to  the  release  point  were  developed  in 
Appendix  C. 

The  first  method  derives  the  control  law  from  the  magnitude  equality  and  the  sec¬ 
ond  method  derives  the  control  law  from  the  direction  equality.  Other  methods  were 
suggested  by  the  AFAL.  One  of  them  involves  using  the  v  and  w  components  of  the 
vectors  instead  of  the  magnitude  and  direction  to  solve  for  the  TOF  and  the  steering 
commands. 

These  control  laws  must  be  evaluated  in  the  control  law  design  phase  and  com¬ 
bined  with  the  FIREFLY  bombing  algorithm  to  design  an  integrated  generalized  mechani¬ 
zation  for  bombing. 


3.  5  MEASUREMENT  FUNCTION  AND  AVIONICS  SUBSYSTEM 


The  measurement  function  receives  as  its  input  the  outputs  of  all  possible  avionic, 
subsystems  xmder  consideration  by  the  AFFA.  Associated  with  each  avionic  sensor 
output  is  a  requirement  to  specify  all  errors  of  quantities  measured  by  that  sensor 
(such  as  means,  variances,  and  power  spectra)  to  the  extent  that  these  are  known 
from  avionic  sensor  test  data.  The  measurement  function  must  also  contain  generalized 
'mbdng"  logic  such  that  simple  mixing  of  data  (data  fusion)  from  sensor  outputs  can  be 
selected  and  supplied  as  inputs  to  the  state  estimator  function  via  a  suitable  choice  of 
the  measurement  specialization  parameter,  Msp.  The  specification  and  design  of  Msp 
will  enable  any  meaningful  combination  of  sensor  quantities  to  be  selected  for  a  given 
avionic  configuration.  Figure  24  shows  the  avionics  subsystem /measurement 
function/generalized  mechanization  interface.  The  measurement  function  supplies  the 
observation  vector  y,  and  the  collection  of  statistical  information  I  such  as  observa¬ 
tion  error  covariance,  correlation  time,  measurement  biases,  etc.  to  the  remainder 
of  the  GM.  In  some  applications  information  must  also  be  supplied  directly  from  the 
measurement  function  to  the  state  estimator  in  order  to  turn  the  Kalman  gains  on  or  off 
In  accordance  with  the  availability  of  measurements.  Ihe  measiurement  function  will 
also  have  to  make  any  coordinate  transformations  necessary  on  sensor  data. 

Table  4  shows  the  avionics  subsystem /measurement  function  (MF)  interface. 
Under  the  avionics  subsystems  (left  column  in  the  table)  all  applicable  avionics  sensors 
are  listed.  The  outputs  of  these  sensors  become  inputs  to  the  MF  which  generates  the 
generic  inputs  to  the  GM  (outputs  of  the  MF).  The  measurement  specialization  param¬ 
eter,  which  must  be  specified  by  the  user  of  the  GM,  determines  what  sensor  combina¬ 
tion  is  available  and  thixs  determines  the  specific  function  to  be  performed  by  the  MF. 
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TABLE  4.  AVIONICS  SUBSYSTEM/MEASUREMENT  FUNCTION  INTERFACE 


AVIONICS  SUBSYSTEM 
(Measurement  Function  Input) 

INPUT  TO  REMAINDER  OF  GM 
(Measurement  Function  Output) 

INS 

Ownshlp  Inertial  Position 

• 

Velocity  (Vj^,  V^,  V^) 

•  Acceleration  (A^.A^.A^) 

•  Glmbal  Angles  (Hdg,  0,  0) 

•  Lat  -  Lon.  (♦.A) 

Ownshlp  Inertial  Velocity  (V^) 

• 

Ownshlp  Inertial  Position 

Ownshlp  Inertial  Acceleration  (Ag) 
(Sensed) 

STRAPDOWN  IMU 

• 

Acceleration  (A  ,  A  ,  A  ) 

U  V  w 

• 

Euler  Angles  (i|/,  9,  0) 

Ownshlp  Transformation 

From  Inertial  to  Body 

• 

Body  Rates  (p,  q,  r) 

Ownshlp  Angular  Rates  u) 

AIR 

DATA  SENSOR 

• 

True  Air  Speed  (V.j.^g) 

True  Airspeed  (V^^g) 

• 

Angle  of  attack  (a) 

Angle  of  Attack  (a) 

• 

Sideslip  Angle  (B) 

Angle  of  Sideslip  (B) 

• 

Relative  Pressure  (p/p  ) 

0 

(if  measurement  is  available) 

Relative  Air  Density  (p/p^) 

Relative  Air  Density  (pl Pq) 

• 

Temperature  (T) 

• 

Pressure  Altitude  (h) 

• 

Mach  Number 

RADAR  (DOPPLER) 

• 

Velocity-ground  speed  (V^) 
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TABLE  4.  AVIONICS  SUBSYSTEM/MEASUREMENT  FUNCTION 
INTERFACE  (CONTINUED) 


AVIONICS  SUBSYSTEM  INPUT  TO  REMAINDER  OF  GM 

(Measurement  Function  Input) _ (Measurement  Function  Output) 

TARGET  RELATIVE  STATE  VECTOR 


ELECTRO-OPTICAL  TRACKER 

•  Range 

•  Glmbal  Angles  (Tg,  Eg) 

•  Tracker  Rates  (01^.,  ui^) 

•  Glmbal  Angle  Rates (Tg,  Eg) 
RADAR  TRACKER 

•  Range  (R) 

•  Glmbal  Angles  (Tg,  Eg) 

•  LOS  Rates  (o).  ,  (o  ,  10.) 

t’  c  d 

•  « 

•  Glmbal  Angle  Rates  (Tg,  Eg) 

•  Range  Rate  (R) 

INS 

•  Glmbal  Angles  (9,  (i,  Kdg) 

BARS 

•  Euler  Angles  (9,  Ql,  Hdg) 

MUZZLE  VELOCITY  SENSOR 

•  Muzzle  Velocity  (Vj^) 

HELMET  SIGHT 

•  Helmet  Angles  (T^,  Ej^) 

or  Direction  Cosines 


•  Range  (R) 

•  Range  Rate  (R) 


•  Tracker  Angular  Rates 

•  Tracker  axis  direction  cosines 

GRAVITY  VECTOR 

•  Transformation 
Platform  to  Body 

e  Transformation  Vertical  Gyro 

Platform  to  aircraft  body. 

GUN  MUZZLE  VELOCITY  (V„) 

M 

HELMET  ORIENTATION 

Direction  Cosines  or 
Helmet  to  Body  Transformation 


RADAR  ALTIMETER 
•  Altitude  above  surface. 


ALTITUDE  (h) 


4.0  GENERALIZED  MECHANIZATION  REQUIREMENTS 

The  primary  objective  of  the  Advanced  FIREFLY  Assessment  (AFFA)  Program  • 
Is  to  provide  a  single  software  tool  that  Is  versatile  enough  to  simulate  specific  mech¬ 
anizations  and  to  evaluate  the  designs  of  these  Integrated  Fire  Flight  Control  (IFFC) 
system  mechanizations  for  a  wide  variety  of  aircraft  and  avionic  subsystems.  In  this 
regard,  design  considerations  for  the  GM  must  Include  the  Interdependency  of  different 
requirements  of  the  functions  that  make  up  the  GM  portion  of  the  AFFA. 

The  gross  structure  of  the  GM  has  been  identified  In  Section  2  and  expanded  In 
detail  In  Section  3  of  the  GMRR.  This  section  identifies  specific  requirements  within 
each  of  the  four  major  functions  comprising  the  GM: 

1.  Fire  Control  Solution  Function 

2.  State  Estimator  Function 

3.  Control  Laws  Function 

4.  Measurement  (Avionic  Subsystem  Inputs)  Function 

In  addition,  requirements  for  coordinate  frames  and  a  definition  of  symbology 
to  be  used  for  the  GM  design  are  Included  In  this  section. 

In  general,  the  GM  requirements  can  be  divided  into  three  categories: 

1.  Functional  Requirements 

2.  Accuracy  Requirements 

3.  Hardware  (or  equipment)  requiren\ents 

The  requirements  category  mostly  dealt  with  In  this  GMRR  are  the  functional 
requirements.  The  accuracy  requirements  will  be  addressed  in  subsequent  tasks  of 
the  AFFA  program.  However,  a  top  level  requirement  is  identified  herein,  that  can 
be  stated  as  follows: 

"The  generalized  mechanization  and  each  of  its  constituent  functions 

shall  meet  the  accuracy  requirements  imposed  upon  it,  such  that  the 

overall  specified  weapon  delivery  accuracy  is  met. " 

4.1  FIRE  CONTROL  FUNCTION  REQUIREMENTS 

The  generalized  mechanization  (GM)  for  Firefly  shall  compute  the  fire  control 
solution  for  three  basic  weapon  delivery  modes:  AAG,  AGG  and  bombing.  The  GM 
shall  be  based  upon  but  not  limited  to  Firefly  II  algorithms  and  shall  encompass  other 
than  Firefly  weapon  delivery  modes  listed  on  Table  5  on  the  following  page. 
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•TABLE  5 

WEAPON  DELIVERY  MODES 


Alr-to-Air 

.  20  mm  gun  director  submode 

Gunnery 

.  LCOS  submodes  for  non- tracked  targets 

(AAG) 

.  Transition  from  non-track  to  track  modes 

Alr-to-Ground 

.  All  CCRP  types  of  weapon  delivery  for 

Gunnery 

(AGG) 

projectiles  and  rockets 

Bombing 

.  Direct  and  Indirect  sub  modes 

.  CCRP/INS  direct  submode 

.  Angle  rate  bombing  system 

.  Indirect  (or  blind)  coordinate  bombing 

.  Offset  bombing 

.  Stick  bombing 

Moreover,  the  GM  of  the  fire  control  solution  algorithms  shall  satisfy  the 
following  additional  requirements; 

1.  Fire  control  algorithms  shall  be  expressed  In  terms  of  the  vectors  used 
In  the  fire  control  geometry.  This  approach  allows  for  modifications  to 
be  made  In  any  of  the  vector  expressions  without  affecting  the  overall  fire 
control  structure,  and  hence  Is  consistent  with  the  GM  requirements. 

2.  Versatility  shall  be  provided  for  using  either  closed-form  expressions  or 
munerlcal  Integration  techniques  for  computing  the  weapon  ballistic 
trajectory. 

3.  The  weapon  TOF  computation  algorithm  shall  converge  rapidly  to  produce 
a  fire  control  solution  In  real  time  consistent  with  airborne  computer  com 
putatlonal  constraints. 

4.  The  effect  of  sight  parallax,  gun  parallax,  bomb  parallax,  and  bomb  rack 
delay  shall  be  Included  In  the  derivation  of  fire  control  equations. 

5.  Fire  control  equations  shall  be  developed  for  the  large  -angle  case  without 
making  the  small  angle  approximation  made  In  the  FIREFLY  II  report. 


6.  The  Generalized  Mechanization  of  ballistic  algorithms  shall  be  capable  of 
controlling  the  delivery  of  the  weapons  listed  below.  The  fire  control  solu¬ 
tion  shall  have  the  caps’  ty  of  calculating  the  Impact  point  and  time  of 
flight  of  all  these  weapons. 

20mm  Guns 
30mm  Guns 
Rockets 

Low  Drag  General  Purpose  Bombs 
Two-stage/ Cluster  Weapons 

7.  Numerical  Integration  algorithms  shall,  to  the  maximum  extent  possible, 
be  Identical  for  the  AAG,  AGG  and  bombing  modes.  Further,  the  Integra¬ 
tion  algorithms  shall  be  highly  adaptable  to  different  types  of  ungulded 
ordnance. 

8.  The  fire  control  solution  shall  compute  the  following  variables: 

Variables  Common  To  All  Weapon  Delivery  Modes 

•  Future  target  position 

•  Predicted  weapon  Impact  point 

•  Predicted  miss  distance 

•  Weapon  TOF  from  release  to  Impact 

•  Ballistic  trajectory  component  vectors  along  the  total  velocity  vector  and 
along  gravity. 

Additional  Variables  for  the  AAG  and  AGG  Modes 

•  Gunllne  direction  cosines  for  a  hit 

•  Desired  range  for  a  hit 

•  Gun  lead  angle 

•  Aircraft  angular  rate  commands. 

Additional  Variables  for  Bombing 

o  Weapon  release  point 

•  Time  to  go  before  release 

•  Steering  commands 

•  Range  vector  R  from  current  attacker  position  to  a  point  directly  above 

P 

the  desired  weapon  Impact  point  on  the  turning  plane. 
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Additional  Variables  for  Stick  Bombing 


•  Distance  to  go  along  a  linear  flight  path  before  initiating  the  turn  and  the 
associated  turning  rate. 

•  Range  vector  R  defined  above  for  each  of  the  endpoints  of  a  linear  target. 

P 

4.2  STATE  ESTIMATION  FUNCTION  REQUIREMENTS 

The  state  estimators  shall  consist  of  ownshlp  state  estimator,  atmospheric 
estimator,  and  target  state  estimator.  In  addition  to  these  estimators  the  state  esti¬ 
mation  function  shall  predict  the  future  target  position  and  the  kinematic  lead  vector 
by  using  the  target  state  estimates. 

4.2.1  Ownshlp  State  Estimator  Requirements 

The  ownshlp  state  estimator  final  requirements  must  necessarily  await  an  error 
budget  derived  from  detailed  error  analysis  of  the  postulated  available  measurements. 
Since  ownshlp  rotational  quantities  such  as  angular  rates  and  orientation  are  usually 
measured  accurately  by  sensors,  there  may  be  no  need  to  estimate  them. 

The  ownshlp  state  estimator  for  the  GM  design  shall  be  based  on  the  following 
assumption  which  becomes  a  requirement  on  the  OFP. 

"The  avionic  subsystem  will  Include  an  INS  or  a  strapdown  system  which 
measures,  at  a  minimum,  the  ownshlp  Inertial  acceleration." 

From  a  qualitative  and  structural  point  of  view  the  following  GM  estimator  re¬ 
quirements  can  be  defined. 

4. 2. 1. 1  Inertial  Acceleration  and  Velocity  Estimator 

The  Inertial  acceleration  and  velocity  estimator  shall  perform  the  following 
functions: 

1.  Estimate  the  ownshlp  inertial  acceleration  and  velocity  from  measurements 
provided  by  glmbaled  and/or  strapdown  IMU's. 

2.  Be  capable  of  computing  ownshlp  inertial  acceleration  and/or  velocity  by 
fusing  any  linear  combination  of  acceleration  or  acceleration  Integrals 
respectively. 
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3.  Compute  Inertial  velocity  from  doppler  radar  estimated  velocity  and  earth 
surface  velocity.  (See  Section  3.3.1. 1) 

4.  Provide  any  linear  combination  of  velocity  estimates  for  Improving  accuracy 
or  testing  degraded  modes. 

All  transformations  required  to  generate  a  useful  ownshlp  motion  output  In  any 
appropriate  coordinate  system  shall  be  defined  and  provided. 

4. 2.1. 2  Angle  of  Attack  and  Sideslip  Estimators 

The  a,  y9  estimator  shall  provide  accurate  estimates  of  air  velocity  vector  direc¬ 
tion  cosines  or  estimates  of  the  angle  of  attack  and  sideslip  angle,  from  pertinent 
measured  or  computed  quantities.  These  quantities  may  Include  the  following; 

1.  Sensed  body  accelerations  and  velocities  from  the  IMU. 

2.  Aircraft  body  rates. 

3.  a  and/or  y9  vane  measurements,  If  any. 

The  o,  P  estimators  shall  be  capable  of  producing  accurate  estimates  for  large 
aircraft  maneuvers  (large  a  close  to  stall  and  large  sideslips).  The  estimates  shall 
not  depend  critically  on  the  aircraft  stability  and  control  derivatives. 

4.2.2  Atmosphere  Estimator  Requirements 

The  atmospheric  estimator  shall  have  the  capability  to: 

1.  Estimate  the  ownshlp  relative  air  velocity, 

2.  Obtain  smooth  estimates  of  air  density,  atmospheric  pressure  and  baro 
altitude  above  mean  sea  level, 

3.  Compute  the  wind  velocity  relative  to  Inertial  coordinates  using  the  esti¬ 
mated  ownshlp  Inertial  velocities  and  the  estimated  air  velocities. 

4. 2. 3  Target  State  Estimator  Requirements 

Essential  to  any  gunfire  control  system  Is  the  ability  to  predice  the  future  tra¬ 
jectory  of  the  target  during  the  TOF  of  the  projectile.  Future  target  position  Is  pre¬ 
dicted  from  current  target  state  estimates  obtained  by  processing  LOS  and  ownership 
measurements  through  an  Extended  Kalman  filter. 
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In  addition  to  providing  accurate  estimation  of  the  target  states,  the  generalized 
target  state  estimator  must  satisfy  the  following  requirements: 

1.  Interface  with  various  avionic  sensors  Including  multiple  sensors. 

2.  Be  compatible  with  various  weapon  delivery  modes,  such  as  AAG,  AGG 
and  bombing. 

3.  Be  compatible  with  alternate  (inertial  and  LOS)  coordinate  frames  for 
estimating  target  states. 

4.  Have  a  modular  design  with  Interchangeable  computational  blocks;  such  as 
alternate  target  acceleration  subfunction  and  alternate  Kalman  gain 
subfunction. 

5.  Be  based  on  a  representative  target  acceleration  model  of  moderate  com¬ 
plexity  which  Includes  acceleration  rate  dependence  on  target  acceleration 
as  well  as  on  target  velocity. 

6.  Have  a  means  of  correcting  the  Kalman  gains  In  the  event  modeling  errors 
tend  to  produce  filter  divergence. 

7.  Account  for  state  dependent  measurement  noise  sources  In  the  filter  design. 

8.  Compensate  for  correlated  measurement  noise  and  measurement  bias, 
whenever  possible.  In  the  estimator  design. 

9.  Any  augmented  state  variables  which  represent  measurement  bias,  drift 
rates,  misalignment,  etc. ,  shall  be  fed  back  to,  or  accounted  for  In  the 
estimator  design  to  Improve  the  state  estimation  performance. 

10.  The  target  state  estimator  shall  be  as  invariant  as  possible  for  various 

tactical  aircraft  with  differing  avionics  and  shall  provide  maximum  common¬ 
ality  for  specialized  mechanizations. 

In  order  to  design  a  practical  and  implementable  target  state  estimator  which 
satisfies  all  of  the  requirements  listed  above,  it  Is  necessary  to  make  the  following 
simplifying  assumptions. 

1.  All  avionic  sensor  outputs  are  synchronized  and  updated  at  the  same  rate 
through  an  "avionic  Interface"  before  they  are  supplied  to  the  generalized 
estimator.  Additional  errors  and  correlation  introduced  by  measurement 
averaging  will  be  accounted  for  as  part  of  the  measurement  function  design. 
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2.  Sensor  output  mixing  (or  data  fusion)  Is  performed  In  the  "measurement 
function"  exter  .al  to  the  state  estimator.  Data  compression  considerations 
and  pertinent  coordinate  transformation  will  be  Included  In  the  measurement 
function  design. 

3.  Signals  required  to  drive  the  target  tracker  will  be  generated  externally 
to  the  state  estimator. 

4.2.4  Future  Tai'get  Position  Prediction  and  Kinematic  Lead 

Target  position  one  TOF  In  the  future  shall  be  computed  using  target  state  esti¬ 
mates,  the  computed  TOF,  and  representative  target  acceleration  models. 

4.3  CONTROL  LAW  REQUIREMENTS 

The  primary  purpose  of  an  aircraft  control  law  is  to  supply  controlled  Inputs 
to  the  flight  control  system  (FCS).  In  conformity  with  the  GE  type  of  control  laws  the 
Inputs  to  the  FCS  shall  be  in  the  form  of  pitch,  yaw,  and  roll  rate  steering  command 
Inputs  and  they  shall  be  compatible  with  various  weapon  delivery  modes  (AAG,  AGG, 
and  bombing)  and  for  differing  conventional  and  CCV-type  aircraft.  Alternate  weapon 
delivery  mode  and  aircraft  compatibility  Is  provided  in  the  FCS  Interface  (or  coupler) 
which  conditions  the  steering  command  Inputs  before  they  are  applied  to  the  FCS. 

The  objectives  of  the  control  law  are  twofold: 

•  To  generate  steering  command  inputs  for  various  weapon  delivery  modes, 
and 

•  To  adapt  the  steering  conunand  inputs,  through  a  coupler,  to  a  specific 
aircraft. 

The  coupler  shall  be  designed  based  on  the  following  requirements: 

1.  The  combination  of  coupler  and  FCS  shall  provide  wideband  and  well 
damped  aircraft  responses  to  angular  rate  steering  command  inputs  as 

the  aircraft  characteristics  for  a  given  flight  condition  would  permit.  This 
requirement  may  be  met  by  compensating  the  FCS  for  a  specific  aircraft 
to  achieve  the  desired  responses. 

2.  The  coupler  shall  be  designed  to  meet  the  gain  and  phase  margin  stability 
requirements  of  the  augmented  aircraft  specified  by  MIL-F-94901. 
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3.  Proportional  plus  Integral  contaroller  shall  be  used  in  the  forward  loop  to 
provide  unity  transfer  function  in  the  steady  state  as  shown  in  Figure  23. 

4.  Authority  sharing  between  the  pilot  and  the  automatic  Firefly  steering 
system  shall  be  built  into  the  coupler  design  with  f.tlot  override  provision. 

5.  The  automatic  steering  command  signals  and  aircraft  load  factor  shall  be 
limited  for  a  specific  aircraft  to  ensure  structural  integrity  and  to  prevent 
the  aircraft  from  getting  into  an  uncontrollable  flight  regime. 

6.  Control  parameters  of  the  coupler  (gains  and  time  constants)  shall  not  de¬ 
pend  directly  on  aircraft  stability  derivatives  and  flight  control  parameters, 
but  on  aircraft  design  limits  such  as  load  factor  and  roll  rate. 

7.  Wherever  possible  the  coupler  shall  be  designed  in  such  a  way  as  to  mini¬ 
mize  the  sensitivity  of  aircraft  responses  to  aircraft  dynamic  parameter 
variations  in  the  flight  envelope. 

8.  The  FCS  shall  have  gust  rejection  capability.  Attitude  sensitivity  to  gust 

is  desirable  in  gunnery  whereas  load  factor  insensitivity  to  gust  is  desirable 
for  bombing. 

9.  Full  utilization  of  the  direct  lift  and  direct  side  force  capability  of  the  CCV- 
t)T)e  aircraft  shall  be  made  in  the  coupler  design  to  improve  the  offensive 
and  survivability  capability  of  the  aircraft. 

10.  Decoupling  control  laws  shall  be  incorporated,  especially  for  the  AAG  mode 
as  indicated  in  Section  3.4.?.  to  eliminate  or  minimize  interaction  between 
longitudinal  and  lateral  aircraft  responses. 

11.  These  additional  control  laws  shall  not  produce  imdesirable  response 
characteristics  a’thou^  some  compromise  in  stability  may  be  inevitable. 

In  many  instances  the  use  of  decoupling  control  laws  and/or  gust  allevia¬ 
tion  control  systems  would  lead  to  neutral  stability. 
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In  the  11  July  1978  correspondence  to  the  AFAL/HWT  it  was  stated  that: 

"Design  efforts  for  the  Control  Augmentation  System  (CAS) 
coupler  are  deleted  as  suggested.  Northrop  will  utilize  the 
F- 15/16  CAS  designs  defined  by  the  Firefly  II  studies  and 
assume  that  CAS  designs  for  the  A- 10,  ATF,  and  CCV  aircraft 
will  be  provided  by  the  Air  Force.  Northrop  will  work  with 
the  Air  Force  to  establish  the  simplified  flight  control-aerodynamic 
transfer  function  representatives  based  upon  the  specified  CAS 
designs. 

If  a  necessity  develops  during  the  course  of  the  study  to  improve 
CAS  characteristics  to  optimize  the  SAM  avoidance  concept,  the 
Air  Force  will  be  informed  and  additional  design  efforts  may  be 
established  at  a  future  time.  " 

Iherefore,  the  coupler  design  for  a  variety  of  aircraft  Is  outside  the  scope  of 
the  present  AFFA  Program.  It  is  assumed  that  the  user  of  the  generalized  Firefly 
mechanization  will  supply  the  appropriate  coupler  in  conjunction  with  the  specific  air¬ 
craft  model  to  be  simulated. 

Requirements  which  relate  to  the  generation  of  steering  commands  which  apply 
dire  ctly  to  the  GM  design  are  listed  below: 

1.  The  steering  control  laws  shall  be  updated  at  such  a  rate  as  to  prevent 
aliasing  and  Interaction  with  the  structural  modes  but  not  faster  than 
necessary  to  conserve  computing  time. 

2.  The  control  law  shall  include  the  proper  compensation  such  as  rate  mixing 
to  compensate  for  the  aircraft  response  lag,  the  computation  time  in  the 
fire  control  computer,  and  the  update  time  on  the  MUX  bus.  Compensation 
shall  not  affect  closed-loop  stability  and  shall  preferably  be  done  outside 
the  FCS  feedback  loops. 

3.  Bombing  control  laws  shall  Include  provision  for  stick  bombing  and  shall 
Include  the  flexibility  of  using  either  the  roll  or  the  pitch  axis  control 
system  of  a  generic  aircraft  to  allow  the  Firefly  control  laws  to  be  Imple¬ 
mented  on  the  A-10  with  minimum  modifications. 

The  control  law  requirements  listed  above  apply  to  fixed  guns  and  fixed  stores  but  not 
to  trainable  guns.  For  fixed  guns,  fuselage  aiming  for  AAG  Is  achieved  by  rotating 
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the  aircraft  in  pitch  and  yaw  relative  to  earth  until  the  total  projectile  velocity  is 
pointed  to  the  desired  aimpolnt  which  tates  into  account  weapon  ballistics  and  gravity 
drop. 

4.4  MEASUREMENT  FUNCTION  AND  AVIONICS  SUBSYSTEM  REQUIREMENTS 

The  distinction  between  the  avioni^-’  subsystem  and  the  measure  ment  function 
is  clearly  defined  in  Section  3. 5  of  this  GMRR.  The  avionic  subsyctem  shall  include 
all  applicable  sensors  and  shall  supply  the  measurements  to  the  measurement  function 
(MF).  The  latter  shall  generate  all  generic  parameters,  such  as  range,  ownshlp 
velocity  and  orientation,  etc. ,  required  for  the  fire  control  solution.  The  inputs  and 
outputs  of  the  MF  shall  include ,  but  not  be  limited  to  the  measurement  quantities  listed 
in  Table  4.  Selection  of  appropriate  combinations  of  sensors  or  measurement  shall 
be  provided  by  the  measurement  specialization  parameter  "Msp".  Specification  of 
Msp  shall  provide  the  required  generic  inputs  to  the  GM  from  a  specific  combination 
of  sensors. 

The  requirements  for  the  MF  are  listed  below: 

1.  The  MF  shall  provide  measurement  mixing  and  data  fusion  as  required 
to  generate  a  single  output  quantity  for  each  generic  variable  such  as 
range,  aircraft  velocity,  etc. 

2.  All  outputs  of  the  MF  shall  be  synchronized  and  updated  at  the  same  rate, 

3.  Additional  measurement  errors  introduced  as  a  result  of  measurement 
averaging  (due  to  correlation)  for  data  synchronization  shall  be  accounted 
for  in  the  measurement  error  covariance  computation. 

4.  Each  sensor  shall  be  identified  as  having  or  not  having  any  internal 
filtering,  smoothing,  or  state  estimation.  For  each  sensor  output  quan¬ 
tity  or  set  of  quantities  which  have  been  filtered,  there  shall  be  specified 
a  measurement  error  covariance  matrix,  bias  error,  and  a  filtering  lag 
time. 

5.  State  dependent  noise  and  correlated  measurement  noise  characteristics 
of  each  sensor  shall  be  identified  and  this  information  shall  be  supplied 
to  the  state  estimators  for  inclusion  in  the  design  process. 

6.  Data  compression  shall  be  accounted  for  and  pertinent  coordinate  trans¬ 
formation  shall  be  provided  in  the  MF. 
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4. 5  COORDINATE  FRAMES  REQUIREMENTS 


The  coordinate  systems  used  for  the  fire  control  solution  shall  be  orthogonal  and 
ri^t-handed  and  shall  consist  of  the  following: 

N,  E ,  D  -  Inertial  or  earth-fixed  coordinates; 

N  -  north,  E  -  east  and  D  -  down, 
u ,  V ,  w  -  Attacker  body  coordinates, 

u  -  forward  along  the  roll  axis,  v  -  out  to  right  wing 
along  the  pitch  axis,  and  w  -  down  along  the  yaw  axis, 
u.  ,  V  ,  w  -  Attacker  wind  coordinates, 

u^  -  along  attacker  velocity  relative  to  the  air  mass. 
v„,  -  to  the  ri^t  and  w„  -  down. 

s,  I,  m-  Sight  coordinates, 

s  -  along  the  target  line-of-sight, 

I  -  out  to  the  right  and  m  -  down. 

t ,  c ,  d ,  -  Tracker  coordinates,  t  -  along  the  tracker  axis, 

c  -  to  the  right  along  the  tracker  elevation  axis  and 
d  -  completing  the  orthognnul  system. 

®h  *  '^h  *  °^h  “  Helmet  si^t  coordinates,  Sj^  -  along  the  target  si^t, 

4  out  to  (he  right  and  down. 

Any  vector  can  be  transformed  from  one  set  of  coordinates  to  another  set  of  coor¬ 
dinates  by  using  the  following  transformation  matrices. 
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Attacker  Body  to  Sight  Transformation  [T] 
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Attacker  Body  to  Helmet  Sl^ht  Transformation  [H] 
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The  generalized  mechanization  shall  use  the  most  efficient  means  of  constructing 
the  transformation  matrices  [E],  [W],  [T],  [F],  (Hj  and  others,  if  any.  Quaternions, 

Euler  angles,  direction  cosines  shall  be  among  the  candidate  systems  to  be  considered 
in  the  design.  The  selection  between  these  alternate  means  of  computing  the  transfor¬ 
mation  matrices  should  be  based  upon  physically  measurable  quantities  available  on  a 
specific  aircraft. 

4.6  NOMENCLATURE  REQUIREMENTS 

1.  The  GM  shall  be  designed  using  the  nomenclature  defined  in  this  section. 

The  coordinate  transformations  and  coordinate  components  shall  be  as  de¬ 
fined  in  Section  4. 5 

2.  The  fire  ccmtrol  mode  point  and  vector  definitions  shall  be  those  described 

in  Figure  25  for  AAG  and  AGG,  Figure  26  for  the  turning  plane  geometry 
for  bombing,  Figure  27  for  the  vector  diagram  at  turn  initiation  for  bombing, 
and  Figure  28  for  the  vector  diagram  at  the  release  point  for  bombing. 

3.  Vectors  shall  be  represented  by  symbols  with  an  over  bar. 

4.  Estimated  parameters  shall  be  represented  by  symbols  with  hat. 

5.  Time  derivatives  will  be  represented  by  overdots,  one  for  first  and  two  for 
second. 
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Vectors 


Parallax,  attacker  IhJU  to  sight 

Parallax,  attacker  IMU  to  gun 

Present  range  to  target 
Kineiiintic  lead  vector 
Gravity  drop 

Weapon  travel  along  total  velocity  vector  V  +  V 

M  G 

Attacker  velocity  at  gun  station 
Muzzle  velocity  along  gunline 


figure  25.  VECTORS  ASSOCIATED  WITH  AAG  AND  AGO  FIRE  CONTROL  MODE 
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CIRCULAR  FLIGHT  PATH 


C  IMU  location  (or  aircraft  CG) 

at  turn  initiation 

R  Bomb  release  point 

P  Point  on  turning  plane 

directly  above  weapon 
Impact  point 

I  Weapon  impact  point 


Sy  Unit  vector  along  V_. 

Rp  Vector  from  point  C  to ‘point  P 

Rp  Weapon  travel  along  V  from 
release  to  point  P. 

G'  Vertical  drop  of  weapon 

Ro  Radius  of  turning  circle 


Angle  between  and  Rp 


FIGURE  26.  -  TURNING  PLANE 

S  R  T 


KL  Kinematic  lead 
Dr 

or  Target  position  offset 


GEOMETRY  FOR  BOMBING 

POINTS 

C  IMU  location  at  turn  initiation 

S  Target  tracking  sensor 

T  Pross  It  target  position 

F*  Future  target  position 

F  Desired  weapon  impact  point 


FIGURE  27.  VECTOR  DIAGRAM  AT  TURN  INITIATION 
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VECTORS 


W  Wind  velocity  relative  to 
Inertial  coordinates 

Pg  Bomb  parallax 

Vg  Ownship  velocity  at  bomb 
location 

Vg  Bomb  ejection  velocity 
relative  to  body 

^CA*  Defined  above. 


« 

FIGURE  28.  VECTOR  DIAGRAM  FOR  BOMB  RELEASE  POINT  COMPUTATION 
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APPENDIX  A 


FIREFLY  n  CONCEPTS 


1.0  CONCEPT  DESCRIPTION 

FIREFLY  n  offers  a  way  to  couple  the  aircraft's  fli^t  and  fire  control  systems 
in  order  to  introduce  refinements  in  weapon  delivery  through  the  use  of  partial  automatic 
control  of  these  tasks.  In  the  case  of  A/A  combat,  the  Integration  of  flight  and  fire 
control  systems  results  in  improvemeits  in  gunnery  accuracy  as  well  as  in  reductions 
of  pilot  work  load.  For  A/G  missions,  the  effects  of  ilREFLY  II  are  a  decrease  in  the 
vulnerability  to  anti-aircraft  artillery  (AAA)  and  surface-to-air  missile  (SAM)  threats 
and  an  increase  in  the  weapon  delivery  accuracy.  These  are  consequences  of  the  fact 
that  FIREFLY  offers  an  alternative  to  the  traditional  technique  of  rolling  out  to  a  wings 
level  attitude  before  releasing  the  weapons.  In  contrast  to  this,  FIREFLY  11  allows 
bombing  during  execution  of  a  turn  in  a  plane  which  is  not  restricted  to  be  horizontal. 

For  AGG,  FIREFLY  facilitates  the  execution  of  an  evasive  maneuver  while  simultaneously 
maintaining  the  gun  line  on  the  target.  The  possibility  of  evasive  aircraft  maneuvers 
while  engaged  in  A/G  combat  increases  the  problem  for  the  AAA  and  SAM  by  making  the 
flight  path  less  predictable  and  the  maneuver  requirements  more  difficult. 

The  Integration  of  fire  and  fli^t  control  systems  is  achieved  in  FIREFLY  II  by  the 
generation  of  command  signals  which  are  used  to  steer  the  aircraft.  The  system  has 
been  designed  with  the  capability  to  offer  several  levels  of  pilot  participation  in  the  con¬ 
trol  tasks,  varying  from  pilot-aided  semi-automatic  control  to  manual  control. 

For  automatic  control  of  AAG,  the  basic  philosophy  in  the  FIREFLY  II  design  is 
to  induce  colinearity  of  the  angular  velocities  of  the  aircraft  and  of  the  LOS  to  the  tar¬ 
get's  future  position  and,  at  the  same  time,  to  zero  out  the  angular  errors  of  the  gun. 

In  AGG  the  basic  approach  in  FIREFLY  EL  is  to  augment  the  sideslip  angle  during 
the  firing  interval  rather  than  to  coordinate  the  aircraft,  as  is  done  traditionally.  The 
presence  of  a  sideslip  angle  does  not  cause  significant  aiming  errors  if  its  value  is  esti¬ 
mated  properly.  It  does  lead  to  a  reduction  in  the  time  available  for  an  AAA  to  achieve 
an  effective  burst,  as  well  as  an  increase  in  the  miss  distance  for  AAA  fire  directed 
£^ainst  the  aircraft. 
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The  FIREFLY  II  bombing  system  offers  the  capability  to  deliver  weapons  in 
high-G  turns  initiated  from  any  dive  or  roll  angle,  and  does  not  call  for  precise  con¬ 
trol  of  aircraft  lift.  The  approach  consists  of  determining  the  release  conditions  as  a 
fimction  of  the  turn  rate  magnitude.  The  range  at  which  the  bomb  should  be  released 
is  displayed  to  the  pilot,  whose  task  is  to  continue  to  turn  at  a  constant  rate  until  the 
release  range  is  reached. 

In  the  following  paragraphs,  various  topics  relating  to  the  General  Electric 
FIREFLY  n  concept  will  be  discussed.  These  topics  Include  the  familiarization  with 
the  HREFLY  n  concepts,  methods  for  determining  the  weapon  delivery  accuracy  when 
FIREFLY  n  is  implemented,  and,  finally,  work  performed  in  the  analysis  of  the 
FIREFLY  n  control  laws  using  the  Terminal  Aerial  Weapon  Delivery  Simulation 
(TAWDS)  program. 

Although  fire  control  algorithms  were  developed  in  the  General  Electric  FIREFLY 
II  report,  considerable  work  was  done  by  Northrop  during  the  familiarization  phase  to 
extract  the  Information  and  tie  the  various  parts  of  the  algorithm  together  into  an  inte¬ 
grated  fire  control  system, 

2.0  FIREFLY  R  FAMILIARIZATION 

The  block  diagram  of  the  FIREFLY  n  system  for  the  AAG  and  ACG  weapon  de¬ 
livery  modes  is  shown  in  Figure  A-1.  Symbols  used  in  this  Figure  and  in  the  equations 
that  will  follow  are  defined  in  Table  A-1. 

As  shown  in  Hgure  A-1,  the  angular  rate  of  the  antenna,  the  measured  target 
range,  and  the  off-boresight  angle  of  the  LOS  to  the  garget,  as  well  as  ownship  velocity 
and  acceleration  are  inputs  to  the  target  state  estimator.  The  estimates  of  the  target 
variables  are  then  used  to  generate  the  target  relative  state  vector,  which  is  required 
as  input  to  the  gunnery  eqxiations.  These  equations  compute  the  desired  angular  rates 
that  should  be  imparted  to  the  gim  and  also  supply  estimates  of  the  gim-aiming  errors. 
The  outputs  of  the  gunnery  equations  serve  to  generate  the  pitch,  yaw,  and  roll  rate 
commands  which,  in  turn,  affect  the  aircraft’s  position,  velocity,  and  acceleration  and 
effectively  close  the  loop. 
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FIGURE  A-1.  FIREFLY  H  FIRE  CONTROL  SYSTEM  FOR  AAG  AND  AGG 


TABLE  A-1.  DEFINITION  OF  TERMS  FOR  FIREFLY  U  CONTROL  SYSTEM 


‘kJG*  ^CG 

®M 

V 


^R 


‘^shn 

e...  e„ 

V  W 

Pi  q,r 


T, 


=  Gun  pitch  and  yaw  rate  commands 
=  Off-boresight  angle  of  target  LOS 
=  Aircraft  velocity 

=  Target  velocity  relative  inertial  coordinates 
=  Relative  target  velocity 
=  Target  range  vector 
=  Range  rate 
=  Measured  target  range 
=  Target  acceleration 
=  Aircraft  body  acceleration 
=  Relative  target  acceleration 
=  Ownshlp  acceleration 
=  Angular  rate  of  tracker  p-^s 
=  Angular  rate  of  target  LOS 
=  Angvilar  velocity  of  tracker  coordinates 
=  Aiming  error  components 
=  Roll,  pitch,  yaw  rate  of  aircraft 
=  Process  noise  in  target  model 
=  Future  target  range 
=  Attacker’s  airspeed  and  muzzle  velocity 
=  Time  of  flight 


a,0  =  Angle  of  attack,  sideslip 


«s!!s<«amf!iassss:n 


The  preceding  description  ondts  the  fact  that  the  coordinate  systems  associated 
with  each  functional  block  are  not  the  same  and  that  transformations  from  one  system 
to  the  other  are  required  throu^out  the  loop.  These  coordinate  systems  include  iner¬ 
tial  coordinates,  aircraft  body  coordinates,  (uvw),  LOS  coordinates  (slm)  and  roll 
stabilized  tracker  coordinates  (scd). 

A2.1  Target  State  Equations 

The  following  differential  equations  define  the  target  state  estimator: 


slm 


(0  -0) 
t  8 


5  [f  =<  -  a)  -  2R 

"t  -  “slm  *  * 


(A-1) 


‘t  ■'  +  T| 

scd  'scd 


slm 


a.  -  0)  ,  X  V 
t  slm  t 


A2.2  Target  Relative  State  Equations 

The  target  relative  velocity  along  the  antenna  line  of  sight  is  given  by; 

Vrs  =  « 

The  relative  velocity  normal  to  the  line  of  sight  is  computed  from 
u  _  R  X  (Vj.  -  V) 


The  target  relative  acceleration  is 

^  - 


-  a. 


{A-2) 


(A-3) 
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A2.3  Gunnery  Equatlona 


The  gunnery  equations  include  calculation  of  the  bullet  TOF,  and  of  the  gun 
rate  and  error  signals  which  serve  to  generate  the  aircraft  commands. 


a.  Time  of  Flight  Computation 

Target  position  T^  seconds  in  the  future  is  predicted  using  the  following 


equation. 


R  +  Vj  Tj  +  i  Tf^ 


(A-4) 


The  relationship  between  T^  and  future  range  is  expressed  by 


'  lv\l 


’'a  ^ 


(A-5) 


where  is  the  attacker's  airspeed,  Vj^  is  the  average  muzzle  velocity  and  Cg  is 
the  ballistic  coefficient  which  is  proportional  to  the  square  root  of  total  projectile 
velocity  (Vj,  +  Vj^).  Manipulation  of  the  above  two  equations  leads  to  a  cubic  poly¬ 
nomial  in  Tj ,  which  is  solved  by  iteration. 

b.  Gxm  Rate  Commands 

The  FIREFLY  control  laws  require  that  the  turn  rate  commands  nominally 
drive  the  gun  at  the  angular  rate  of  the  future  position  of  the  target. 

In  vector  form,  the  gun  angular  rate  commands  Wqq  are  given  by 

[  U  X  +  Tj  +  T^)  +  (U  X  I)  +  T^T^CU  x  7^)  j 


u, 


CG 


where 


-  U  X  (oj  X  AU) 


T  V 
f  M 


(A-6) 


1  4- 


1  + 
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and  where  II  is  a  unit  vector  along  the  aircraft  u-bo(fy  axis,  u  is  the  aircraft  angular 
velocity,  dU  is  the  vector  differences  between  U  and  a  unit  vector  along  the  gunline, 
and  Tj  is  the  time  rate  of  change  of  given  by: 

•  _  R  Tf  C^) _ 

'  2  Tj  (Cj  +0.5  a^g)  +  Cj  R  5’  R  -  -  Vjj  8ln^W/2) 


c.  Required  Range  for  a  Hit 


The  angular  gun  errors  are  defined  as  the  angular  separation  between  the 
LOS  to  the  target  and  the  LOS  required  for  a  hit,  Rj^.  The  latter,  expressed  in  vector 
form,  is  given  by 


1  + 


Vm  - 


1  + 


\  -  I  V  (A-8) 


where  all  symbols  have  been  identified  in  the  preceecting  paragraphs. 

If  both  Rj^  and  the  actual  range  to  the  target  R  are  expressed  in  the  attacker’s 
body  axes  (u,  v,  w),  the  components  about  the  v  and  w  axes  of  the  gun  pointing  errors 
become: 


R_  .  -  R 
Rw  w 

'lV  ■  - R - 

V  -  "v 

LW  R 

d.  Generation  of  Aircraft  Commands 

Figure  A-2  summarizes  the  main  features  of  the  AAG  control  system. 

Now  q^Q  and  r^Q  are  the  bofiy  axes  components  of  oJqq,  the  angular  rate 
of  predicted  future  target  position,  and  (p,  q,  r)  represent  the  ownship  bo<fy  rates. 

The  rate  commands  are  generated  as  a  linear  weighted  sum  of  the  angular 
velocity  of  the  future  position  of  the  target  and  the  present  gun  aiming  errors. 
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HGURE  A-2.  AIR-TO-AIR  GUNNERY  CONTROL  SYSTEM 
3.0  FIREFLY  II  TAWDS  ANALYSIS 

Following  the  initial  familiarization  with  the  FIREFLY  II  control  laws,  a  simu¬ 
lation  of  their  performance  was  carried  out  using  a  six-degree  cf  freedom  nonlinear 
program  named  Terminal  Aerial  Weapon  Delivery  Simulation  (TAWDS),  (Reference  6), 
Both  A/A  and  A/G  versions  of  this  program  exist,  but  only  the  A/A  FIREFLY  n  laws 
were  Implemented. 

An  abbreviated  block  diagram  of  TAWDS  (A/A)  Is  shown  In  Figure  A-3,  As  Indi¬ 
cated,  the  program  simulates  an  attack  on  a  target  by  a  weapon  system  consisting  of  a 
pilot  and  a  fighter  aircraft  with  Its  fll^t  control,  sight,  and  we^on  delivery  systems. 
The  TAV/DS  program  propagates  stationary  and  dynamic  source  errors  Into  statistical 
Impact  error  alstrlbutlons.  Stationary  source  errors  considered  arise  from  the  flight 
profile  and  system  mechanization  errors.  Dynamic  source  errors  are  those  due  to 
atmospheric  disturbances,  weapon  release  forces,  and  pilot  steering  tasks.  The  per¬ 
formance  output  data  for  the  TAWDS  (A/A)  program  are  time  history  responses  for 
the  attacker  and  target  aircraft  responses,  plots  of  specified  attack  aircraft  we^on 
delivery  time  history  responses,  stationary  source  error  printouts,  ensemble  pass 
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statistics  printouts  for  the  tracking  error  responses,  and  ensemble  pass  statistics 
printouts  for  the  bullet  Impact  error. 

The  approach  taken  to  Incorporate  the  FEREFIY II  control  laws  Into  TAWDS  was 
to  replace  Ihe  blocks  Inside  the  dashed  rectangle  in  Figure  A-3  equivalent  ones,  as  . 
shown  In  Figure  A-4.  By  this  choice,  the  aircraft  was  flown  exclusively  using  the 
commands  generated  by  FIREFLY  II  rather  than  be  a  mbcture  of  pilot  and  automatic 
control  commands.  Substitution  of  the  aircraft's  control  system  dynamics  and  airframe 
response  by  simplified  transfer  functions  were  dictated  by  the  fact  that  the  aircraft 
modeled  In  TAWDS  Is  an  F-4  and  that  the  detailed  aerodynamic  data  required  by  TAWDS 
was  not  available  for  the  F-15  aircraft.  For  this  reason,  transfer  functions  relating 
the  rate  commands  to  the  actual  bocty  rates  v/ere  derived  using  data  from  the  F-lS's 
FIREFLY  n  results.  An  additional  simplification  was  introduced  by  assuming  that  the 
velocity,  angle  of  attach  and  sideslip  angle  of  the  attacker  were  constant  during  the 
encounter.  As  was  the  case  in  the  FIREFLY  II  report,  a  director  gunslght  was 
considered. 

The  results  of  the  simulation  showed  that  the  elevation  and  traverse  tracking 
errors  stayed  below  0. 5  degrees  In  an  encounter  between  an  F-15  and  a  maneuvering 
target  represented  by  an  F-4. 
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FIGURE  A-3.  TAWDS  (A/A)  BLCXJK  DIAGRAM 
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FIGURE  A-4.  FIREFLY-TAWDS  (A/A)  INTERFACE 


APPENDIX  B 


DATA  FUSION  TECHNIQUES 


Techniques  for  fusing  multiple  sensor  outputs  are  summarized  and  the  relation¬ 
ship  between  these  techniques  and  the  Advanced  FIREFLY  Assessment  program  are 
indicated  in  Table  B-1. 


TABLE  B-1 


RELATIONSHIP  OF  DATA  FUSION  TECHNIQUES  TO  AFFA  PROGRAM 


Data  Fusing  Techniques 

Relationship  to  AFFA  Program 

Selection  of  one  sensor  out 
of  multiple  sensor 

Integral  Part  of  FIREFLY  II  GM 

• 

Optimally  mixed  sensors  based  on 
-  Noise  Covariance 

Noise  spectrum 

Time  response 

New  Concepts 

• 

Seml-adaptlve  measurement 
mixing 

• 

Parallel  Filtering 

• 

Filtered  Sensor  Output 

• 

Fully  adaptive  Kalman 
filtering 

Outside  the  scope  of  AFFA  Program 

The  simplest  means  of  dealing  with  multiple  sensors  which  measure  the  same 
variable  is  to  pick  the  more  accurate  of  the  two  for  the  particular  operating  condition. 
This  scheme  was  considered  as  part  of  Advanced  FIREFLY  GM.  On  the  other  hand,  a 
fully  adaptive  Kalman  filter  estimator  was  considered  to  be  outside  the  scope  of  the 
AFFA  program.  The  remaining  data  fusion  techniques  listed  in  Table  B-1  were  classi¬ 
fied  as  "new  concepts"  and  were  investigated  during  the  requirements  analysis  phase 
of  the  AFFA  Program. 


Optimally  mixing  the  sensor  outputs  improves  the  accuracy  of  the  measurement 
(or  the  estimate),  with  minimum  addition  of  the  fire-flight  control  system.  For  in¬ 
stance,  these  techniques  may  provide  between  80  percent  and  90  percent  of  the  accuracy 
achieved  through  a  fully  adaptive  Kalman  filter,  at  a  fraction  of  the  complexily.  More¬ 
over,  the  fully  adaptive  Kalman  filter  usually  has  a  slower  response  time  than  an 
ordinary  Kalnaan  filter  to  prevent  instability  caused  by  interaction  between  the  filtering 
loop  and  the  loop  which  adaptively  varies  the  Kalman  gains  in  response  to  the  measure¬ 
ment  error  covariance.  The  flow  filter  time  respons*^  is  an  undesirable  feature  of  the 
fully  adaptive  Kalman  filter. 

Combining  multiple  sensor  outputs  provides  back  up  protection  in  the  event  of 
sensor  failure  and  sensor  jamming  in  addition  to  reducing  the  measurement  error. 

Since  the  radar  is  easily  jammable,  a  jamming  detector  cai'  be  used  to  detect  the 
presence  of  jamming  signals  and  switch  to  another  sensor  automatically.  Despite  the 
advantages  of  data  fusion  there  are  limited  applications  wherein  the  available  sensor 
outputs  can  be  combined  optimally.  Topical  examples  include  mixing  the  angle  measure 
ments  from  EO  sensor  and  monopulse  radar  and  combining  tho  range  measurements 
from  radar  and  the  laser. 

Among  the  techniques  listed  above  in  Table  B-1  those  techniques  which  can  be 
easily  included  in  the  GM  Involve  adding  multiple  sensor  outputs  with  fixed  gains. 

Other  techniques  are  sensor-dependent  ana  therefore,  they  do  not  conform  with  the 
GM  concepts.  These  more  complicated  techniques  could  be  studied  in  follow-on  work 
aa  they  relate  to  specific  applications  and  specific  combinations  of  sensor. 

Optimal  Mixing  Based  on  Measurement  Err  or  Covariance 

Figure  B-1  shows  a  technique  for  mixing  tho  multiple  sensor  outputs  based  on 
apriori  knowledge  of  measurement  noise  covariances.  The  fixed  gains,  and  Cg  in 
this  figure  are  given  by 

R, 

’  ^2  *  R^  +  R^  (B-l) 

where 

Rj^  ■  Cov(Xj^)  ;  R^  =  CovCX^) 
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FIGURE  B-1.  OPTIMAL  MIXING  BASED  ON  MEASUREMENT  ERROR  COVARIANCE 
Hie  error  covariance  of  the  optimally  mixed  output,  X, 

\  <=-"> 

is  smaller  than  the  covariance  of  either  sensor  output.  Hence,  by  mitring  the  sensor 
outputs  as  shown  in  Figure  B-1  the  measurement  error  covariance  will  be  reduced 
and  the  accuracy  of  the  measurement  will  be  Improved.  Maximum  improvement  is 
obtained  v/hen  Rj^  =  Rg  in  which  case  the  covariance  of  X  is  one  half  the  covariance  of 
either  sensor  output  as  shown  by  the  knee  in  the  curve  in  Figure  B-2.  As  the  ratio 
Rf/Rg  deviates  from  unity,  the  improvement  decreases  to  the  point  where  sensor  out¬ 
put  mixing  no  longer  becomes  practical.  For  instance,  to  decrease  the  covariance  to 
80  percent  of  the  error  covariance  of  the  more  accurate  sensor,  ^1/^2 
tween  0.25  and  4.  CXitside  this  interval  mixing  the  sensor  outputs  becomes  impractical. 


If  the  measurement  error  covariances  Rj^  and  Rg  are  not  known  precisely,  setting 
the  gain  to  other  than  result  in  suboptimal  mixing  in  which  the 

combined  measurement  is  still  more  accurate  than  either  measurement,  as  long  as 
lies  in  the  interval: 


1 


R2  +  R1 


(B-3) 


for  Rj^  <  Rg*  Outside  this  interval,  the  mixed  output,  X,  will  be  less  accurate  than  the 
more  accurate  sensor  output  X^^.  This,  of  course,  defeats  the  purpose  of  data  fusion. 
Therefore,  ballpark  values  for  the  meiisurement  error  covariances  must  be  known 
before  attempting  data  fusion. 
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SENSOR  1 


SENSOR  2 


APPENDIX  C 


A-10  CONTROL  LAW  MODIFICATION 


1.  INTRODUCTION 

FIREFLY  II  algorithms  generate  a  roll  control  lav;  for  automatically  steering 
the  aircraft  toward  the  release  point.  The  roll  axis  control  system  on  the  A-10  Is  not 
accessible  for  Inserting  external  commands  without  major  modifications  to  the  flight 
control  system.  Therefore,  the  feasibility  of  using  the  pitch  axis  control  system  of 
the  A-10  (which  is  accessible)  Instead  of  the  roll  axis  control  system,  for  automati¬ 
cally  steering  the  aircraft  toward  the  release  point  is  investigated  herein  as  one  of 
Uie  "new  concepts"  In  the  AFFA  program. 

Incorporation  of  a  modified  control  law  for  the  A-10  into  the  FIREFLY  II  bomb¬ 
ing  algorithm  will  simply  Involve  switching  the  appropriate  control  axis  and  using  the 
the  associated  computational  method. 

2.  FIREFLY  R  BOMBING  CONCEPT 

The  FIREFLY  11  bombing  system  provides  for  a  maneuvering  approach  to  bomb 
release  from  a  non-wxngs  level  attitude.  There  exists  a  fundamental  relationship 
among  the  encounter  parameters  to  execute  the  proper  maneuver  for  accurate  bomb 
release.  This  relationship  takes  the  form  of  equating  the  magnitude  and  direction  of 
two  vectors  associated  with  the  bomb  release  geometry. 

The  basic  guidance  equations  for  bombing  are  obtained  from  Figure  C-1  which 
shows  the  geometry  on  a  turning  plane.  From  this  geometry  one  has  the  following 
fundamental  relationship  tliat  must  exist  among  the  encounter  parameters  if  the 
attacker  is  to  execute  a  constant  rate  maneuver  from  present  position  to  the  release 
point  at  a  specific  TOF: 


'R 


V,  R 

a  p 


sin  X 
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(C-l) 


Present 

Attacker 

Position 


- . ..A 


Release 
•  Point 
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(Center) 


T 

Target 


Point  directly  above  target 

Center  of  circular  bombing  trajectory 

Radius  of  circle 

Attacker  velocity  vector  relative  to  air  mass 
Angular  rate  vector 
Unit  vector  along  V 
Unit  vector  along  R 


FIGURE  C-1.  TURNll^G  PLANE  GEOMETRY  FOR  FIREFLY  BOMBING 


Using  vector  notation  Equation  (C-1)  ii  n'^.^orltten  as  follows: 


S  X  R 

V  p 


2  2 
%  -  V 


(C-2) 


Equation  C-2  represents  the  vector  solution  to  the  correct  release  condition  and 
the  vectors  on  both  sides  of  the  equation  must  be  equal  in  magnitude  and  direction. 

The  approach  used  in  the  FIREFLY  II  mechanization  is  to  treat  the  magnitude  of  the 
aircraft  turn  as  an  independent  variable  and  find  the  bomb  TOF  (and  the  release  point) 
which  makes  the  magnitudes  of  the  two  vectors  C  and  a»  equal.  The  direction  of  the 
two  vectors  Cp  and  <j  is  then  made  equal  by  rotating  the  angular  rate  vector  by  control¬ 
ling  roll  rate,  which  steers  the  aircraft  toward  the  release  point. 
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Ordinarily,  the  pilot  controls  the  aircraft  turning  rate  by  moving  the  stick 

laterally  until  the  desired  bi:nk  angle  (0)  or  the  desired  yaw  rate  (r  =  sln^ )  Is 

achieved,  after  which  the  stick  is  rebirned  to  the  neutral  position.  Therefore,  In  a 

steady  turn,  the  roll  rate  about  the  wind  axis  Is  zero,  and  hence  In  wind  coordinates 

the  vectors  S  ,  w  and  C  are  defined  by 

V  p  , 


1 

'0  ' 

'Si' 

'o  ■ 
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;  oi  = 

S? 

S2 
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0 
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w_ 

_S3_ 

_S2_ 

The  pilot  can  also  command  a  pitch  rate  by  moving  the  stick  longitudinally  If  he 
so  desires  to  make  a  bombing  approach  on  a  nonhorizontal  turning  plane.  The  vector 
sum  of  pitch  rate  (q^)  and  yaw  rate  (r^)  determines  the  magnitude  of  the  turning 
rate  {u).  The  TOF  is  computed  such  that  the  magnitudes  of  the  vectors  in  Equation 
(C-2)  are  made  equal;  and  that  is: 

R  ^  2 

ISl  - 

The  roll  control  law  is  obtained  from  the  requirement  that  the  directions  of  the 
vectors  and  w  be  made  equal.  Since  both  vectors  and  w  are  in  the  plane  perpendi¬ 
cular  to  the  wind  axis,  the  angular  error  between  those  vectors  is 


e 


r 

w 


(C-5) 


The  error  <  vanishes  when  the  vectors  C  and  w  are  colirear.  Therefore,  the 

P 

FIREFLY  n  control  strategy  is  to  feed  the  error  <  into  the  roll  axis  control  system  In 
such  a  way  as  to  null  the  error  e  as  shown  In  Figure  C-2  below. 


Because  of  the  Inherent  Integration  between  the  commanded  roll  rate  and  the 
achieved  bank  angle,  the  error  €  will  be  nulled  and  the  direction  of  the  vectors  C  and  w 
v;ill  be  made  equal  in  the  steady  state,  by  controlling  yaw  rate  (r^)  via  aircraft  body 
roll  rate  (p). 


3.  PROPOSED  ROLL/PrrCH  CONTROL  LAW  MODIFICATIONS  FOR  THE  A-10 

Tiere  are  two  possible  approaches  for  adapting  the  FIREFLY  n  bombing  control 
laws  to  the  A-10.  The  first  approach  retains  the  basic  FIREFLY  II  control  philosophy 
while  using  the  pitch  axis  control  system  (instead  of  the  roll-axis  control  system)  to 
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FIGURE  C-2.  FIREFLY  H  ROLL  AXIS  CONTROL  SYSTEM 

Steer  the  aircraft.  The  second  approach  reverses  the  functions  of  steering  control 
and  TOF  computation  as  described  below. 

First  Approach 

•  The  pilot  sets  the  magnitude  of  the  turn  rate  |cJ1  by  specifying  a  yaw  rate. 

•  TOF  Is  computed  such  that  the  magnitudes  of  the  vectors  and  T5  are 
made  equal. 

•  The  direction  of  w  is  adjusted  automatically  by  controlling  pitch. 

The  equations  for  computing  the  TOF  unmodified.  The  aircraft  control  law, 
however,  Is  modified  as  shown  In  Figure  C-3„ 

Notice  that  a  proportional  plus  Integral  controller  Is  needed  In  the  pitch  control 
loop  In  order  to  achieve  essentially  zero  error  €  and  ensure  that  the  directions  of  the 
vectors  C  and  w  are  the  same  In  the  steady  state. 

Second  Approach 

•  In  this  approach  the  pilot  still  sets  the  approximate  magnitude  of  |u| 
by  adjusting  yaw  rate. 
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FIGUBE  C-3.  PROPOSED  PITCH  AXIS  CONTROL  SYP’.i  NO,  1  FOR  THE  A-10 

•  Pitch  axis  control  system  adjusts  the  magnitude  of  u  until  the  magnitudes 
of  the  vectors  in  Equation  (C-4)  are  made  equal. 

•  TOF  is  then  computed  such  that  the  direction  of  the  two  vectors  C^ 
and  Zi  is  made  equal. 


From  Equation  (C-4)  one  hasj 


l“l  ■  +  <l„^  +  r  ^ 


i'v’-'S.l 


w 


2  2 


2  V 


(C-6) 


Since  the  ri^ht  hand  side  is  independent  of  fee  magnitude  of  T5,  equality  can  be 
achieved  by  adjusting  q^.  The  pitch  axis  control  system  for  the  second  approach 
takes  the  form  of  Figure  C-4. 

Because  of  the  square  root  operation,  fee  control  system  in  Figure  C-4  should 
2  2  2 

not  be  engaged  unless  Z  >  p  +  r  •  This  condition  can  be  achieved  by  simply 

W  t/ 

decreasing  the  turn  rate  magnitude. 

The  relationship  for  computing  the  TOF  is  then  obtained  i?y  nulling  the  error  £ 
in  Equation  (C-5).  A  Newton-Raphson  algorithm  similar  to  the  one  used  in  the 
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FIGURE  0-4.  PROPOSED  PITCH  AXIS  CONTROI,  SYSTEM  NO,  2  FOR  THE  A-10 

FIREFLY  II  mechanization  can  be  used  for  this  purpose.  The  recursive  equation 
for  TOF  Iss  (e) 


at. 

(C-7a) 

L  fJ 

n 

e  = 

^2  * 

r 

w 

(C-7b) 

3e  '^2 

•  Cl  — 

^3 

-  -  -  •  t: 

(C-7c) 

1 

rt 

at,  w 

f 

The  partial  derivatives  aR^g/atj  ,  eRpg/at^  must  be  evaluated  in  closed-form  or 
numerically  before  a  solution  can  be  obtained  for  t^.  The  convergence  properties  of 
this  recursive  equation  must  be  verified  through  simulation. 


4.  FUTURE  WORK 

Other  methods  for  adapting  the  FIREFLY  II  bombing  algorithms  to  the  A-10 
should  also  be  investigated,  and  the  most  promising  approach  should  be  selected  on 
the  basis  of  a  simulation  study.  One  method  which  appears  feasible  is  to  equate  the 

components  of  vectors  and  U.  y — YVk — }  ^w  instead 


of  the  direction  and  magnitudes  of  the  respective  vectors. 
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